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ABSTRACT 


Heat  transfer  and  hydrodynami c  performance  of  three  dif- 
ferent spirally  fluted  tubes  was  determined.   The  tubes  were 
5/8"  in  nominal  diameter  and  were  made  of  aluminum.   Results 
were  compared  to  5/8"  OD,  smooth  copper-nickel  and  aluminum 
tubes . 

Data  was  taken  by  condensing  steam  at  about  3  psia  on  the 
outside  surface  of  a  hori zontal ly  mounted  tube  in  the  center 
of  a  tube  bank.   The  center  tube  was  cooled  by  water  on  the 
inside  at  velocities  of  3  to  25  feet  per  second.   The  overall 
heat  transfer  coefficient  was  determined  directly  from  exper- 
imental data.   The  inside  and  outside  heat  transfer  coeffi- 
cients were  determined  using  the  Wilson  plot  technique.   The 
cooling  water  pressure  drop  was  measured  inside  the  tube  and 
converted  to  the  friction  factor  in  the  enhanced  section. 

The  overall  heat  transfer  coefficients  of  the  enhanced 
tubes  were  as  large  as  1.75  times  the  corresponding  smooth 
tube  yalue  for  the  same  mass  flow  rate  of  cooling  water.   The 
inside  heat  transfer  coefficients  increased  by  about  a  factor 
of  3  while  the  outside  heat  transfer  coefficients  decreased 
by  10  to  20  percent  when  compared  to  smooth  tube  values. 

The  results  of  this  work  indicate  that  the  required  con- 
denser surface  area  can  be  reduced  by  50  percent  if  these  en- 
hanced tubes  are  used  in  place  of  smooth  tubes. 
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I.   INTRODUCTION 

A.   BACKGROUND  INFORMATION 

In  recent  years,  there  has  been  an  increased  awareness 
regarding  the  use  of  enhanced  heat  transfer  surfaces  in  the 
design  of  heat  exchangers.   By  using  enhanced  heat  transfer 
surfaces,  heat  exchangers  can  be  designed  to  be  smaller  in 
size  which  can  result  in  a  savings  in  capital  costs.   In 
addition,  they  can  also  be  designed. to  be  more  efficient 
which  will  save  on  operating  costs. 

A  recent  study  by  Search  [1]  of  an  air  aircraft  carrier 
main  condenser  concluded  that  marine  condenser  heat  loads,  at 
constant  pumping  power,  can  be  increased  by  up  to  50  percent 
using  heat  transfer  enhancement  techniques.   He  also  concluded 
that  the  use  of  corrugated  copper-nickel  tubes  can  be  used  to 
decrease  condenser  weight  and  cost  for  constant  pumping  power 
and  constant  heat  load  when  compared  to  the  conventional  con- 
denser design.   Search  also  concluded  that  size  and  weight 
savings  on  the  order  of  40  percent  could  be  realized  depending 
on  the  heat  transfer  enhancement  method  used. 

The  possibilities  of  such  savings  has  spurred  much  research 
and  many  design  efforts  in  recent  years.   Bergles  [2,3]  has 
summarized  extensive  works  in  both  single  phase  and  two  phase 
heat  transfer  enhancement.   He  has  compiled  and  listed  the 
many  efforts  that  have  been  undertaken.   In  addition,  Bergles 
has  compared  the  results  of  several  different  experimenters 
using  a  variety  of  performance  criteria. 
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In  the  particular  case  of  steam  condensers,  a  variety  of 
tests  have  been  conducted  to  study  the  behavior  and  performance 
of  enhanced  tubes.   Palen,  Chan,  and  Taborek   [4]  tested 
several  Turbotec  tubes  manufactured  by  Spiral  Tubing  Corpor- 
ation.  The  purpose  of  these  tests  were  to  compare  the  per- 
formance of  the  Turbotec  tube  to  the  smooth  tube.   The  tests 
were  conducted  using  a  tube  bundle  consisting  of  196  tubes. 
Steam  at  55  psig  or  105  psig  was  on  the  shell  side  of  the 
condenser,  while  water  at  mass  flow  rates  between  380,000  Ibm/hr 
and  800,000  lbm/hr  were  used  on  the  tubeside.   The  average 
bulk  temperature  of  the  cooling  water  was  245  to  320  °F. 
EissenbeYg  L5]  performed  an  extensive  study  of  condenser  tube 
heat  transfer  coefficients  using  a  multi-tube  bundle.   He 
conducted  his  tests  on  copper-nickel  roped  tubes  using  a  steam 
vapor  temperature  between  160  and  300   F.   Newson  and  Hodgson 
[6]  conducted  heat  transfer  experiments  on  32  different  types 
of  tubes.   All  of  their  experiments  were  conducted  with  satur- 
ated steam  at  atmospheric  pressure  (212°F)  on  the  outside  of 
a  single  vertically  mounted  tube.   Water  was  pumped  through 
the  inside  of  the  tube  at  a  rate  sufficient  to  maintain  the 
LMTD  constant  at  11°F.   Watkinson  et  al .  [7]  conducted  tests 
on  18  Noranda  Forge  Fin  tubes.   Their  tests  were  conducted 
with  a  single  tube  orientated  in  a  horizontal  position.   Steam 
was  supplied  to  the  outside  of  the  tube  at  atmospheric  pres- 
sure or  at  a  slightly  higher  pressure  while  water  was  used  on 
the  inside.   Catchpole  and  Drew  [8]  conducted  experiments  on 
five  radially  grooved  tubes.   In  these  tests  steam  was  supplied 
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at  2  psia  and  the  cooling  water  velocity  was  maintained  at 
10  ft/sec.   The  tubes  could  either  be  tested  in  a  single  tube 
arrangement  or  as  a  bundle.   Young  et  al .  [9]  compared  the 
Korodense  tube,  manufactured  by  the  Wolverine  Division  of 
Universal  Oil  Products,  to  a  smooth  tube.   These  tests  were 
conducted  at  two  different  steam  temperatures  of  100°F  and 
212°F.   The  cooling  water  velocity  was  varied  from  about 
3  ft/sec  to  about  6.5  ft/sec.   Rothfus    [10]  tested  inter- 
nally finned  tubes  made  from  aluminum  and  copper.   The  tubes 
were  mounted  horizontally  and  measured  for  water  side  perform- 
ance only.   The  water  was  heated  by  means  of  heating  tapes 
wrapped  around  the  outer  walls  of  the  test  tube. 

The  tubes  used  in  each  of  the  above  tests  are  widely  dif- 
ferent in  geometrical  shape.   Still  another  tube  design  has 
been  proposed  recently  by  General  Atomic  Company  [11].   It  is 
fabricated  differently  than  the  above  mentioned  tubes,  result- 
ing perhaps  in  different  heat  transfer  performance  when  com- 
pared to  other  enhanced  surface  tubes  already  on  the  market. 
In  addition,  for  each  of  the  experiments  mentioned  above, 
there  was  a  different  test  procedure  and/or  set  of  conditions 
used.   The  need  exists  therefore  for  a  systematic  and  consist- 
ent method  to  test  and  analyze  all  enhanced  surface  tubes. 

Beck  [12]  designed  a  test  facility  at  the  Naval  Postgrad- 
uate School  that  permits  the  testing  of  a  single,  horizontally 
mounted,  condenser  tube.   Pence  [13]  built  and  tested  this 
system.   He  conducted  his  tests  using  a  smooth  copper-nickel 
tube.   The  results  of  Pence's  tests  indicate  that  the 
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test  facility  is  technically  sound  and  that  further  experi- 
ments with  enhanced-surface  tubes  are  possible. 

B.   GOALS  OF  THIS  WORK 

As  mentioned  earlier,  the  enhanced-surface  tube  manu- 
factured by  General  Atomic  Company  was  ready  to  be  tested. 
The  purpose  of  this  thesis  was  then  two-fold: 

1.  to  determine  the  heat  transfer  and  pressure  drop 
performance  characteristics  of  three  different  spirally 
fluted  tubes  manufactured  by  General  Atomic  Company  and  to 
compare  their  performance  to  smooth  tube  operation,  and 

2.  to  develop  a  consistent  method  for  arriving  at  per- 
formance characteristics.   To  accomplish  this,  two  alterna- 
tive data  reduction  schemes  were  used.   One  scheme  was  based 
on  the  smooth  end  geometry  of  the  tube,  which  seemed  to  be 
what  was  used  most  often  in  the  existing  literature.   A  second 
method  was  based  on  the  enhanced  section's  geometry.   To 
accomplish  these  reduction  schemes,  a  computer  program  and  the 
necessary  equations  were  developed. 
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II.   EXPERIMENTAL  FACILITY 

A.  TEST  FACILITY 

The  test  facility  is  seen  in  Figure  1.   The  layout  was 
designed  by  Beck  [12]  and  built  and  tested  by  Pence  [13].   A 
detailed  description  of  the  components  used  in  the  various 
systems  may  be  found  in  these  reports.   Only  a  general  descrip- 
tion of  the  various  systems  will  be  found  within  this  report. 
Particular  attention  will  be  focused  on  the  experimental  tubes, 
their  construction  and  their  location  within  the  test  section. 
See  Appendix  A  for  calibration  procedures  of  components  re- 
quiring calibration. 

B.  STEAM  SYSTEM 

The  steam  system  is  shown  in  Figure  2.   The  boiler  is  an 
electrically  heated  Fulton  Boiler  which  produces  saturated 
steam  at  100  Ibm/hr.   The  steam  leaves  the  boiler  via  a  3/4 
inch  line  and  the  boiler-isolation  valve  (MS-1).   The  water 
contained  in  the  steam  is  removed  by  the  steam  separator.   The 
steam  continues  through  the  system  past  a  flow  meter  and 
through  the  throttle  valve  (MS-3)  where  the  pressure  is  reduced 
The  steam  next  passes  through  the  desuperheater  wherein  water 
from  the  feed  system  is  injected  in  order  to  remove  some  of 
the  sensible  heat  from  the  steam.   The  steam  continues  on  into 
the  test  condenser  where  part  of  it  is  condensed  on  the  test 
tube.   The  steam  not  condensed  is  collected  in  the  vapor  outlet 
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and  sent  to  the  secondary  condenser  wherein  the  latent  heat 
of  vaporization  is  removed.   If  the  boiler  fails,  steam  may  be 
provided  via  the  house-steam-cross-connect  valve  (MS-2).  Steam 
could  be  routed  around  the  test  condenser  to  the  secondary 
condenser  via  the  bypass  valve  (MS-4).   All  steam  lines  (except 
section  downstream  of  MS-3,  see  [13])  were  insulated  with  one 
inch  fiberglass  insulation. 

C.   TEST  CONDENSER 

The  test  condenser  is  shown  in  Figures  3,  4,  and  5.   Steam 
enters  via  the  top.   It  then  passes  through  the  expansion 
section  over  the  baffle  separators,  and  through  three  layers 
of  150  mesh  screen  and  flow  straightener  into  the  tube  bundle. 
The  condensate  collects  at  the  bottom  of  the  test  condenser 
where  it  flows  through  two  1/2  inch  lines  to  the  hotwell. 

The  viewing  window,  shown  in  Figure  3  and  Figure  4,  allows 
viewing  of  the  condensation  process.  Two  types  of  pyrex  glass 
windows  were  used.  One  type  was  a  standard  pyrex  glass  plate, 
1/2  inch  thick.  The  second  type  of  window  was  an  Owens-Corning 
pyrex  glass  with  a  transparent  electrically  conducting  coating 
applied  over  the  surface.  A  Lambda  power  supply  set  at  20  VDC 
and  1.5  A  was  used  with  this  window  to  supply  power  to  the 
electrically  conducting  surface  to  minimize  fogging  effects. 

The  tube  sheet  arrangement  is  as  seen  in  Figure  5.   There 
are  eight  5/8  inch  0D,  18  gauge,  90-10  copper-nickel  tubes 
arranged  in  a  typical  condenser  configuration,  with  a  spacing 
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to  diameter  ratio  (S/D)  of  1.5,  around  a  single  test  tube. 

The  test  tube  is  the  only  tube  with  water  passing  through  it. 
This  arrangement  was  selected  to  best  simulate  the  steam  flow 
conditions  in  an  actual  condenser. 

The  test  condenser  is  insulated  with  two  inches  of  Johns- 
Manville  Aerotube  sheet  insulation. 

D.   CONDENSATE  AND  FEEDWATER  SYSTEMS 

The  condensate  and  feedwater  systems  are  shown  in  Figure 
6.   The  test  condenser  hotwell  collects  the  condensate  from 
the  test  tube,  while  the  secondary  condenser  hotwell  collects 
the  condensate  from  the  secondary  condenser.   Valve  C-l  allows 
isolation  of  the  test  condenser  hotwell  so  the  condensate  mass 
from  the  test  condenser  may  be  measured.   Valve  C-4  is  a  vent 
valve  between  the  test  condenser  hotwell  and  the  test  condenser 
The  condensate  is  pumped  from  the  hotwells  to  the  feedwater 
tank  by  the  condensate  pump.   The  feed  pump  routes  the  water 
from  the  feed  tank  to  the  boiler  via  the  solenoid-controlled 
valve  FW-3,  a  hot-water  filter,  and  the  boiler-isolation 
valve,  FW-4. 

The  feedwater  temperature  is  maintained  between  130  F  and 
140°F  by  thermostat  controlled  heaters.   This  reduces  fluctua- 
tions in  the  boiler  output  and  provides  a  source  of  water  at 
a  temperature  near  saturation  for  the  desuperheater . 

If  house  steam  is  used,  the  condensate  is  returned  to  the 
house  system  via  C-3. 
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The  condensate  lines  are  insulated  with  3/4  inch  Johns- 
Manville  Aerotube  insulation.   The  feedwater  lines  are  insu- 
lated with  1/2  inch  thick  fiberglass  insulation. 

E.   COOLING  WATER  SYSTEM 

The  cooling  water  system  is  a  partially-closed  system  as 
shown  in  Figure  7.   The  water  is  pumped  from  the  supply  tank 
via  a  7-1/2  HP  pump.   The  water  is  routed  to  the  test  tube 
via  one  of  two  rotameters.   A  low  flow  rotameter  allows  up  to 
1.48  gpm  to  flow  through  the  test  tube,  while  a  high  flow 
rotameter  permits  up  to  18.8  gpm.   The  water  returns  to  the 
supply  tank  via  a  dry  cooling  tower.   The  dry  cooling  tower 
was  constructed  using  four  large  radiators  connected  in  series 
The  water  was  directed  through  the  radiators  and  outside  air 
was  forced  over  the  cooling  surface  by  a  centrifugal  fan. 

The  bypass-rotameter ,  downstream  of  CW-3,  is  provided  to 
permit  an  increased  volume  of  water  to  flow  through  the  cool- 
ing tower. 

The  system  piping  was  reduced  from  1  inch  to  5/8  inch 
(approximate  OD  of  all  test  tubes)  at  a  distance  of  approxi- 
mately 2-1/2  feet  ahead  of  the  test  condenser  to  insure  fully 
developed  flow  at  the  test-tube  entrance.   Pressure  taps  were 
installed  in  the  permanent  piping  at  the  ends  of  the  test  tube 
(see  Figure  8)  to  permit  the  measurement  of  the  overall 
pressure  drop. 

The  cooling  water  lines  were  insulated  with  1  inch  thick 
Johns-Manville  Areotube  insulation. 
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F.   SECONDARY  SYSTEMS 

1  .   Vacuum  System 

The  vacuum  is  maintained  by  a  mechanical -vacuum  pump 
and  a  vacuum  regulator  which  induces  an  air  leak  into  the 
vacuum  line.   A  nitrogen  trap  is  provided  at  the  inlet  to  the 
vacuum  pump  to  remove  entrained  water  from  the  vacuum  line 
thus  preventing  contamination  of  the  vacuum  pump  oil. 

2.  Desuperheater 

The  desuperheater  removes  sensible  heat  from  the 
superheated  steam  by  injecting  feed  water  at  about  150°F.   The 
feedwater  flow  into  the  desuperheater  is  controlled  by  DS-1 
and  measured  by  a  rotameter.   The  excess  water  is  collected  in 
a  tank,  located  below  the  desuperheater,  and  returned  to  the 
feed-water  tank  periodically  during  the  experimental  runs. 

3.  Pressure  Tap  (PT)  Tube  Pressure  Drop  Measurement  System 

The  PT  tube  pressure  drop  measurement  system  is  as 
shown  in  Figure  9.   This  system  was  designed  to  allow  connect- 
ing of  the  PT-tube  to  the  manometer  and  to  facilitate  the 
measuring  of  the  various  pressure  drops  across  the  tube.   Figure 
10  shows  how  a  typical  PT  tube  is  connected  in  the  system. 
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G.   INSTRUMENTATION 

1 .  Flow  Rates 

Fulton  rotameters  were  used  to  measure  the  flow  rate 
of  water  in  the  cooling  water  system  and  the  desuperheater , 
while  an  Ellison  Annubar  and  a  differential  water  manometer 
were  used  to  determine  steam  flow. 

2.  Pressure 

Several  different  types  of  pressure  measurement 
devices  were  used  in  this  facility.   They  were:  a  Bourdon 
tube  pressure  gage  which  was  used  to  measure  boiler  pressure, 
a  compound  gage  which  was  used  to  measure  the  secondary  con- 
denser pressure,  an  absolute  pressure  transducer  and  a  30 
inch  mercury  manometer  which  were  used  to  measure  the  test 
condenser  pressure,  and  a  12-foot  mercury  manometer  which  was 
used  to  measure  the  cooling  water  pressure  drop  across  the 
test  tube. 

3.  Temperature 

There  were  three  types  of  thermocouples  used  in  this 
facility.   Stainless  steel  sheathed,  copper-constantan  thermo- 
couples were  used  as  the  primary  temperature  monitoring  devices 
Eleven  temperatures  required  for  data  reduction  were  measured 
using  these  devices.   Table  1  lists  the  locations  monitored. 
Figure  2  shows  the  location  of  six  vapor  space  thermocouples. 
Cooling  water  thermocouples  were  located  as  shown  in  Figure  7. 
Teflon  coated  copper-constantan  thermocouples  were  used  as 
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secondary  measuring  devices.   Table  2  lists  the  locations 
monitored  using  these  thermocouples  and  an  i ron-constantan 
thermocouple  was  used  to  measure  the  boiler  temperature. 

4.   Data  Collection  and  Display 

An  Autodata  collection  system  was  utilized  to  record 
and  display  the  temperatures  in  degrees  Celsius  obtained  from 
the  primary  thermocouples  and  to  record  and  display  the  pres- 
sure in  cm  Hg  inside  the  test  condenser.   See  Table  1  for 
channel  numbers  of  the  temperature  monitoring  devices. 

A  28  channel -di gi tal  pyrometer  was  utilized  to  dis- 
play the  temperatures  obtained  from  the  secondary  thermocouples 
and  a  single  channel  pyrometer  displayed  the  temperature  from 
the  i ron-constantan  thermocouple.   See  Table  2  for  channel 
numbers . 

H.   TEST  TUBES 

The  enhanced  tubes  tested  during  this  study  were  manufac- 
tured by  General  Atomic  Co.   Three  aluminum  tubes  were  made. 
Each  has  an  overall  length  of  48  inches  and  an  enhanced  sec- 
tion length  of  36  inches.   They  have  helical  flutes  on  both 
the  inside  and  outside  surfaces,  which  are  formed  by  running 
a  flat  strip  through  rollers  which  cause  the  flat  surface  to 
become  wavy.   The  wavy  strip  is  then  spirally  wound  and  seam 
welded  to  form  a  tube.   Figure  11  shows  a  short  section  of 
spirally  fluted  tube  rolled  from  strip  and  the  start  of  a 
roll  from  strip.   Figure  12  shows  the  three  tubes  and  their 
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respective  helix  angles  (HA)  manufactured  for  this  experiment. 
Figure  13  shows  a  cross  section  view  of  each  of  these  tubes. 

As  seen  in  Figure  12,  the  helical  portion  of  the  test 
tube,  henceforth  referred  to  as  the  enhanced  section,  is 
welded  to  a  smooth  piece  of  aluminum  tube,  hereafter  referred 
to  as  the  smooth  end.   Each  end  of  the  tube  has  this  smooth 
end  to  facilitate  installation  into  the  condenser  tube  sheet. 

In  addition  to  the  above  tubes,  three  special  tubes  were 
manufactured  by  General  Atomic  Co.  to  permit  the  measurement 
of  the  pressure  drop  throughout  the  test  section.   These 
tubes  are  shown  in  Figure  14.   Three  static  pressure  holes 
were  drilled  in  each  tube's  enhanced  section  by  the  EDM  process 
The  holes  located  at  each  end  of  the  enhanced  section  were 
inspected  with  a  boroscope  by  General  Atomic  Co.  and  found  to 
be  free  of  burrs.   Figure  15  shows  the  relative  location  of 
these  taps  as  well  as  the  static  pressure  taps  in  the  smooth 
ends . 
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III.   EXPERIMENTAL  PROCEDURES 

A.   INSTALLATION  AND  OPERATING  PROCEDURES 
1.   Preparation  of  Condenser  Tubes 

Prior  to  any  run,  the  condenser  tubes  had  to  be  pro- 
perly prepared  to  insure  filmwise  condensation.   The  prepara- 
tion of  the  tubes  was  standardized  in  an  effort  to  minimize 
errors  due  to  different  surface  conditions.   The  wall  thermo- 
couple also  had  to  be  prepared  and  installed  in  such  a  manner 
as  to  reduce  the  possibility  of  introducing  errors. 

The  tubes  manufactured  by  General  Atomic  required  a 
special  installation  procedure  into  the  test  condenser  since 
they  have  a  larger  diameter  along  the  enhanced  section  than 
along  their  smooth  ends.   These  procedures  are  described  in 
detail  in  Appendix  B.   Basically,  the  test  tubes  were  prepared 
by  first  installing  a  thermocouple  to  measure  the  wall  temper- 
ature.  Secondly,  their  exterior  and  interior  surfaces  were 
cleaned  to  insure  proper  wetting  characteristics  and  to  insure 
that  all  deposits  were  removed. 

The  copper-nickel  tube,  tube  4  in  Table  3,  was  prepared 
in  accordance  with  the  procedure  given  in  Pence  [13],  whereas 
the  smooth  aluminum  tube,  tube  5  in  Table  3,  was  prepared  in 
accordance  with  the  procedures  listed  in  Appendix  B. 
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2.   System  Operation  and  Steady  State  Conditions 

Pence  [13]  developed  a  detailed  set  of  operating  pro- 
cedures for  this  system.   They  are  i  ncl  uded,  wi  th  minor  changes, 
in  this  report  as  Appendix  C. 

In  general  it  takes  about  three  hours  from  initial 
light-  off  until  steady  state  conditions  are  established.   The 
feedwater  is  heated  up  to  140  F  by  energizing  the  feed  tank 
heaters  and  recirculating  the  water.   After  installation  of 
the  test  tube  is  complete,  the  vacuum  system  can  be  activated. 
The  data  collection  system  is  programmed,  including  setting 
the  date  and  time  in  accordance  with  reference  [14].   The 
cooling  water  system  is  placed  in  operation.   The  18.8  gpm 
rotameter  is  set  at  about  50  percent  flow  to  allow  adequate 
venting  of  both  legs  of  the  12-foot  manometer.   The  rotameter 
is  then  reset  to  the  lowest  flow  point  for  'system  operation. 
The  steam  system  can  now  be  placed  into  operation. 

Steady  state  conditions  must  be  established  prior  to 
data  collection.   To  determine  this,  two  parameters  were 
monitored.   They  were  the  cooling  water  inlet  temperature  and 
the  steam  vapor  temperature.   The  cooling  water  inlet  tempera- 
ture did  not  rise  more  than  l°F/hr.   The  steam  vapor  tempera- 
ture did  not  vary  more  than  6°F  between  the  six  vapor  thermo- 
couples in  the  condenser  nor  did  the  change  in  temperature  at 
an  individual  thermocouple  exceed  0.5°F.   The  steaming  condi- 
tions and  cooling  water  flow  conditions  remained  constant  while 
establishing  steady  state  conditions. 
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The  time  for  the  system  to  stabilize  was  generally 
about  one  hour  which  is  somewhat  less  than  that  reported  by 
Pence  [13].   The  probable  reason  for  this  is  that  with  the 
addition  of  a  much  higher  capacity  cooling  tower,  the  cooling 
water  temperature  stabilized  much  faster  than  in  his  tests. 
This  is  evidenced  by  the  fact  that  the  cooling  water  inlet 
temperatures  were  on  the  order  of  10  F  lower  for  this  experi- 
ment than  for  those  experiments  performed  by  Pence.   The 
ambient  temperatures,  however,  were  nearly  the  same. 

The  procedure  outlined  above  is  for  a  full  data  col- 
lection run.   There  were  two  other  types  of  runs  conducted. 
Isothermal  pressure  runs  were  conducted  on  the  test  tubes  to 
allow  comparison  of  the  non-isothermal  and  isothermal  friction 
factors.   Also  the  special  pressure  tap  tubes  were  tested  by 
connecting  these  tubes  to  the  PT  tube  pressure  drop  measurement 
system  to  allow  measurement  of  pressure  drops  throughout  the 
enhanced  sections.   Both  types  of  pressure  runs  required  only 
the  cooling  water  system  and  instrumentation  systems  to  be 
energized. 

3.   Maintenance  Procedures 

Periodically  the  systems  required  various  forms  of 
maintenance.   Following  each  run,  the  boiler  received  two 
bottom  blow  downs  to  remove  any  sediment  that  may  have  settled. 
The  supply  tank  in  the  cooling  water  system  required  occasional 
refilling.   Water  treated  by  routing  tap  water  through  a  com- 
mercially rented  resin  bed  was  used  to  reduce  the  amount  of 
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contaminants  in  the  water  that  could  deposit  on  the  tubes. 
The  filter  in  the  feedwater  line  required  changing  approxi- 
mately every    three  months  to  prevent  low  boiler  levels  due 
to  lack  of  feedwater.   The  condenser  glass  window  required 
cleaning  after  approximately  five  runs.   This  was  true  whether 
using  the  heated  glass  or  the  standard  pyrex  glass.   Prior  to 
reinstalling  the  pyrex  glass,  a  light  coating  of  Glycerol 
Reagent  ACS  H0CH2CH0HCH20H,  was  applied  to  the  inside  surface 

to  enhance  the  viewing  of  the  condensation  process. 

B.   DATA  REDUCTION  PROCEDURES 

In  evaluating  the  data  obtained  from  the  heat  transfer 
runs,  two  objectives  were  established.   The  first  of  these 
was  to  present  the  data  in  such  a  way  as  to  make  it  immed- 
iately useful  to  the  designer.   The  second  objective  was  to 
establish  a  reduction  scheme  that  would  allow  the  comparison 
of  enhanced  tubes  based  on  their  actual  internal  surface  areas 

1.   Reduction  Based  on  the  Smooth  End  Diameter,  D. 

To  meet  the  condenser  designer's  needs,  it  was 
felt  that  the  data  should  be  reduced  using  the  smooth  end 
diameter.   This  would  allow  a  direct  substitution  of  an  en- 
hanced tube  for  a  smooth  tube  and  is  especially  important  when 
considering  the  comparison  of  a  wide  variety  of  tube  types. 
In  addition,  a  nominal  area  was  defined.   The  nominal  area  was 
based  on  the  outside  surface  area  of  a  5/8-inch  OD  smooth  tube, 
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Appendix  D,  the  sample  calculations,  is  a  complete 
listing  of  the  equations  used  to  evaluate  the  data.   Appendix  E 
is  a  derivation  of  the  probable  error  in  the  data  reduction 
equations,  followed  by  a  sample  error  analysis  for  tube  2, 
45°HA. 

a.   Overall  Heat  Transfer  Coefficient 

The  method  employed  to  arrive  at  the  overall  heat 
transfer  coefficient  is  straightforward  and  similar  to  that 
employed  by  many  researchers  in  the  past. 

The  heat  transfer  rate  to  the  cooling  water  is 
given  by 


m  cp  (TcQ  -  Tc.) 


(1) 


The  heat  transfer  rate  can  also  be  found  from  the  -overall  heat 
transfer  coefficient  by 


U   An  LMTD 
n  n 


(2) 


where 


LMTD 


(Tv  -  Tci>  -  <Tv  -  Tco> 
T   -  Tc 


(3) 


After  combining  equations  (1),  (2),  and  (3)  it  is  found  that 


m  c 
U   =  E 


n 


n 


T   -  Tc 
in  (^ M 


(4) 
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A  schematic  illustration  of  the  procedures  to  arrive  at  U   is 
shown  in  Figure  16. 

To  remove  the  effect  of  the  tube  wall  material,  a 
corrected  heat  transfer  coefficient  is  found  from 


U 


1 


(5) 


w 


where  R,  is  the  calculated  wall  resistance 
w 


b.   Inside  Heat  Transfer  Coefficient 

The  Nusselt  number  on  the  inside  is  found  from 
the  Sieder  Tate  relationship,  found  in  Holman  [15]  as: 


Nu 


h.D. 


„   D  0.8  DJ/3,  .   ,0.14 
ci  Re    Pr    (^/%) 


(6) 


In  the  above  equation,  c.  is  referred  to  as  the  Sieder  Tate 
constant.   The  remainder  of  the  right  hand  side  of  the  above 

equation  (Re  '   Pr  '  (y/y  )    )  will  be  referred  to  as  the 

w 

Sieder  Tate  parameter,  and  the  procedure  for  arriving  at  this 

value  is  illustrated  schematically  in  Figure  17.   The  Wilson 

plot  is  used  to  arrive  at  the  value  of  the  Sieder  Tate  constant 

The  Wilson  plot  was  developed  in  1915  by  Wilson  [16],  and  has 

been  modified  by  several  researchers  since.   The  procedure  used 

in  this  research  was  developed  by  Briggs  and  Young  [17]. 

The  Wilson  plot  is  merely  a  plot  of  1/Un  versus 

the  inverse  of  the  Sieder  Tate  parameter  which  should  be  a 
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straight  line  when  varying  the  cooling  water  velocity.  The 
reasoning  behind  the  Wilson  plot  can  be  seen  in  the  follow- 
ing development. 

The  overall  heat  transfer  coefficient  can  be 
written  as; 

1 


U 


D  •  , 

Dihi     w    hQ 


(7) 


The  inverse  of  equation  (7)  is 


_1_  =   9_     R   J_ 

Un      Dihi     Rw   hQ 


(8) 


If  (R  +  r—  )  i s  assumed  to  be  constant  and  equation  (6)  is 
w   nQ 

solved  for  h.  in  terms  of  the  Sieder  Tate  parameter,  equation 
(8)  can  be  rewritten  as: 


Do  D  -0.8  D  -1/3  ,  ,  ,-0.14  D 
—t—  Re  Pr  (u/uw)  +  B 
c.k  w 


(9) 


where   B 


w   h 
o 


The  form  of  equation  (9)  is  then  exactly  that  of  a  straight 
line, 

Y   =   MX  +  B    ,  (10) 


where 


Y   =  r±- 


(10a) 
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1 

Sieder  Tate  parameter 


,  and 


cikb 


(10b) 
(10c) 


The  values  of  1/tl   and  the  Sieder  Tate  parameter 
f  n  r 

are  obtained  by  varying  the  water  velocity  and  holding  the 
other  parameters,  such  as  water  temperatures,  steam  vapor 
temperatures  and  condenser  tube  wall  temperature,  nearly  con- 
stant.  When  l/ll   is  plotted  versus  Re"0,8  Pr"1/3  (y/u  )~0-14 

n  w      , 

a  linear  regression  subroutine  [18]  fits  these  points  to  a 
straight  line  and  then  solves  for  the  slope,  M,  and  the  inter- 
cept, B.   Knowing  the  slope,  M,  the  Sieder  Tate  constant,  c . , 
can  be  found  from  equation  (10c).   The  inside  heat  transfer 
coefficient,  h.,  is  then  found  from  equation  (6). 

Once  the  inside  heat  transfer  coefficient,  h  .  , 
is  known,  then  the  Nusselt  number  can  be  solved  for  in 
equation  (6),  to  find  the  Stanton  number, 


St 


Nu 


RePr 


cpg 


(in 


The  cooling  water  properties  (p,  y,  k,  c  . 

Pr)  are  solved  for  as  shown  in  Appendix  D.   Appendix  D  also 

demonstrates  the  procedure  for  arriving  at  the  water  viscosity 
evaluated  at  the  condenser  tube  wall,  y  . 


35 


c.  Outside  Heat  Transfer  Coefficient 

The  outside  heat  transfer  coefficient,  h  ,  can 

o 

now  be  found  from  equation  (7).   Figure  18  schematically 
illustrates  the  various  steps  outlined  above. 

d.  Friction  Factor 


from 


The  friction  factor  for  the  test  tube  is  found 


TS 


(Pb)(APTS)(144)(2gc) 


4(LTS/Dt)  G' 


(12) 


The  variables  used  in  equation  (12)  are  solved 
for  as  shown  in  Appendix  D,  while  the  geometric  constants  are 
found  in  Table  3.   APyc*  however,  is  reduced  from  the  measured 
pressure  drop  outlined  below. 

A  P  T  <-  is  the  pressure  drop  in  the  enhanced  section 

of  the  test  tube.   The  measured  pressure  drop,  AP   as  seen  in 

r  r    m 

Figure  8,  is  taken  over  the  entire  tube  length,  a  distance 
of  51-1/8  inches.   Since  the  enhanced  section  is  only  36 
inches  long,  the  pressure  drop  over  each  of  the  smooth  ends 
must  be  subtracted  off  of  the  measured  pressure  drop.   This 
is  done  by  calculating  the  friction  factor  in  the  smooth  ends 
using: 

f   = 


0'*7?5  for  Re  <  30,000 


or 


Re 
0.046 


Re 


0.2 


for  Re  >  30,000 


(13) 
(14) 
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The  smooth-end-section  pressure  drops  can  then  be  calcualted 
from, 


(fs)(4)(LTS/D.)(G£) 
(Pb)(144)(2gc) 


(15) 


As  can  be  seen  in  Figure  13,  the  cross  sectional 
flow  area  of  the  enhanced  section  of  the  test  tube  is  less 
than  the  cross  sectional  flow  area  of  the  smooth  end  of  the 
tube.   Therefore,  the  water  undergoes  an  expansion  and  a  con- 
traction at  the  exit  and  entrance  to  the  enhanced  section  of 
the  tube.   Associated  with  the  expansion  and  contraction  pro- 
cesses are  certain  irreversible  losses  which  cause  additional 
pressure  drops  to  occur.   These  pressure  drops  must  also  be 
subtracted  off  of  the  measured  pressure  drop  and  are  calcu- 
lated as  shown  in  reference  [19]: 


AP 


e/cn 


PVTS  <Ke  +  Kcn> 


(16) 


Since  the  variations  in  the  contraction  and  expansion  coef- 
ficients K   and  K  are  slight  over  the  range  of  Reynolds 
en      e        3  3 

numbers  used,  an  average  of  these  values  was  used  in  equation 
(16). 


Therefore,  APTS  is  found  using  equations  (15) 


and  (16) 


AP 


TS 


APm  -  AP„  -  AP0. 
m     s     e/cn 


(17) 


and  the  friction  factor  for  the  test  section  is  solved  for 
in  equation  (12). 
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e.   Performance  Criteria 

To  compare  the  enhanced,  or  augmented  tubes  with 
the  smooth  tube,  it  was  necessary  to  use  some  meaningful 
performance  criteria. 

Using  the  Colburn  Analogy,  as  found  in  reference 
[20],  provided  one  such  criterion.   Using  this  analogy,  the 
heat  transfer  performance  is  compared  to  the  friction-factor 
performance  as  seen  in  Appendix  D. 

Bergles  [3]  outlines  several  performance  criteria 
based  on  the  inside  heat  transfer  coefficients  by  solving  for 
the  ratio  of  augmented  to  smooth  tube  areas  while  holding 
various  parameters  constant.   One  such  ratio  is  defined  by 


Nu  /Pr 


0.4 


lu  /Pr 
a 


0.4 


(18) 


which  assumes  that  Q,  m,  D .. ,  T.  and  LMTD  are  constant,  and 

R   .  =  R   +  1/h   =  0.   In  equation  (18)  the  augmented  heat 

transfer  coefficient  h   is  the  value  h.  referred  to  earlier. 

a  l 

In  this  situation  the  flow  velocities  for  the  smooth  and 
augmented  tubes  are  the  same. 

The  area  ratio  defined  by  equation  (18)  does  not, 
however,  take  into  account  the  increases  in  pressure  drop 
and  hence  the  increases  in  friction  factor  caused  by  enhance- 
ment techniques.   The  increase  in  pressure  drop  can  be  included 
when  evaluating  the  performance  of  an  enhanced  tube  compared 
to  that  of  a  smooth  tube.   Bergles  [3]  shows  this  by  defining 
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an  area  ratio  for  constant  pumping  power  as  well  as  for  con- 
ditions as  defined  for  in  equation  (18). 

The  pumping  power  is  given  by: 


(p  v  I  D2)  4  f  cfc]   u2-o 


D'  x2g 


(19) 


{£-)     (rrDL)  fv 
yc 


where  ttDL  is  the  inside  surface  area  for  the  tube  in  question 
By  setting  the  pumping  power  of  a  smooth  tube  equal  to  the 
pumping  power  of  an  enhanced  tube,  it  is  found  that: 


Res  fs 
Rea  fa 


(20) 


Notice,  that  in  this  situation  of  constant  pumping  power,  the 
flow  velocities  and  hence  Reynolds  numbers  are  different  for 
the  smooth  and  the  augmented  tube. 

The  heat  flow  rate  is  given  by: 


hi  A.LMTD. 


(21) 


Since  the  heat  flow  is  also  assumed  to  be  constant  in  both 
the  enhanced  and  smooth  tubes,  the  area  ratio  can  be  found 


A 
a 

A7 


Nu  /Pr 


0.4 


Nu  /Pr 
a 


WA 


(22) 
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Equation  (20)  can  now  be  set  equal  to  equation 


(22)  to  show 


Aa     Nus/Pr0-4     Re3s    f 


Aa     Nu  /Pr0"4 


Rea  fa 


(23) 


If  Nu   is  replaced  in  the  above  equation  by  the 
Dittus  Boelter  relationship  as  found  in  L15]: 


Nu 


0.023  Re  °-8  Pr0'4 
s 


(24) 


and  f   is  replaced  by  equation  (14): 


.046 


Re 


0.2 


(14) 


equation  (23)  can  be  solved  for  the  smooth  tube  Reynolds 
number  in  terms  of  the  augmented  conditions; 


Re 


=\SZ 


2Nu  /Pr 
a 


0.4 


(25) 


In  these  expressions, 
Re.   = 


PD1  vs 


(26) 


Knowing  Re  ,  equation  (23)  can  be  solved  for  the  resulting 
area  ratio. 

Since  a  sizeable  portion  of  the  overall  resistance 
in  a  naval  condenser  could  be  caused  by  the  wall  resistance 
and  the  outside  thermal  resistance,  the  area  ratios  as  defined 
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by  Bergles  [3]  should  be  expanded  to  include  these  external 
resistances.   If  the  heat  flow  is  written  by  equation  (2): 


U   A   LMTD 
n   n 


(2) 


and  a  thin  wall  thickness  is  assumed,  then  the  external  re- 
sistance effects  on  the  area  ratio  can  be  included  in  the 
analysis.   The  wall  thickness  must  be  assumed  to  be  small 
since  the  nominal  area  is  based  on  an  outside  diameter  of 
5/8-inch. 

Invoking  all  of  the  assumptions  made  earlier, 
then  the  results  of  the  constant  pumping  power  case  are  again 


In  addition,  constant  heat  flow  results  in 


(20) 


U 


(27) 


As  before,  these  two  area  ratios  can  be  set  equal,  and  it  is 
found  that; 


Aa     Us     vs  fs 


u 


<h 


(28) 


As  mentioned  by  Search  [1],  for  smooth  tubes,  it  is  found  in 
general  that: 


U 


C  Vv 


(29) 
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where  C  is  a  constant  for  a  given  tube  size,  material  and 
water  inlet  temperature.  Also,  from  equation  (14),  it  is 
known  that 


046 


Re 


0.2 


(14) 


By  choosing  v  .  I)  ,  and  f  and  substituting  equation  (26)  for 
a    a       a 

Re   into  equation  (14),  equation  (28)  can  be  solved  for  the 

corresponding  smooth  tube  velocity,  as: 

3 

^     f 

u 


1 

f   v"  C     pD.  ,  /K    273" 
~  [,   w        (    )      1 


L(Ua)(.046) 


(30) 


Knowing  v  ,  then  equation  (29)  is  solved  for  U   and  equation 
(28)  is  solved  for  the  area  ratio. 

In  selecting  the  values  of  the  constants  to  sub- 
stitute into  equation  (30), the  following  procedures  were 
utilized: 

(1)  U  was  corrected  to  70  F  using  the  procedure 

a 

defined  in  reference  [21], 

(2)  C  was  determined  by  using  the  values  of  U 
and  v  for  the  smooth  copper-nickel  tube  in  run  4,  and  solving 
for  C  in  equation  (29).   The  value  for  C  is  not  a  constant 
over  the  range  of  flows  observed,  therefore,  an  average  value 
of  C  =   251,  was  obtained  and  used.   The  resulting  error  in 

U   at  any  velocity  was  less  than  7  percent. 

(3)  The  dynamic  viscosity  used  was  an  average 
value  of  the  viscosities  at  each  flow  point  used  in  the  area 
ratio  analysis. 
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See  Appendix  F  for  a  summary  of  the  area 
ratio  procedures  and  the  results  for  both  cases  where 

Rext  "  °  and  Rext  '  °' 

2.   Reduction  Based  on  the  Hydraulic  Diameter,  D. 


The  reduction  procedures  for  this  method  were  similar 
to  the  procedures  used  for  the  reduction  based  on  D..   The 
major  obstacle  in  obtaining  meaningful  results  was  in  deter- 
mining the  enhanced  section's  geometry,  and  once  determined, 
how  best  to  apply  it  to  the  available  equations. 

a.   Enhanced  Section  Geometry 

It  is  well-known  how  critical  the  measurement  of 
the  hydraulic  diameter  is  to  heat  transfer  and  fluid  flow. 
There  were  two  alternative  methods  available  and  both  were 
examined.   The  volumes  of  the  enhanced  tubes  were  measured 
and,  by  using  the  length  of  the  tube,  the  cross-sectional 
area  was  obtained.   Knowing  the  cross-sectional  area  and  using 
the  perimeter  information  provided  by  General  Atomic,  a 
hydraulic  diameter  could  be  found.   This  diameter  was  discarded 
in  favor  of  the  diameter  found  by  using  a  planimeter  to  measure 
the  wetted  perimeter  from  an  enlarged  view  of  Figure  13.   The 
hydraulic  diameter  is  then  found  from: 

4A_ 


D, 


Pw 


(31) 
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Similar  problems  were  encountered  in  determining 
the  wall  thickness  and  subsequent  wall  resistance.   See 
Appendix  G  for  details  of  this  procedure. 

b.   Heat  Transfer  Coefficients 


To  introduce  this  geometry  into  the  equations  used 
to  solve  for  the  heat  transfer  coefficients,  it  is  first  ne- 
cessary to  recall  that  the  resistance  to  heat  flow  across  a 
tube  is  equal  to  the  sum  of  the  individual  resistances  as 
shown  earlier.   Therefore, 


1 
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A  U 
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met 
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Aihi 


(32) 


where 


A    =   2-rr  Pw 
o  o 


A.   =   2tt  Pw .  ,   and 


Amet=  ^   Pbar 


(33) 
(34) 
(35) 


By  setting  1/A  U   =  1/A  U  ,  and  multiplying  through  by  A 

allows  us  to  obtain  the  overall  heat  transfer  coefficient 
based  on  the  same  nominal  geometry  that  was  used  in  the 
plain-end  reduction.   The  equation  takes  the  form: 
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A  R 
o  w 


A  U 
n  n 


Aihi 


(36) 


o      met 

Equation  (36)  is  solved  exactly  as  was  equation  (8)  in  section 
B.  l.a.  of  this  chapter  to  obtain  the  Wilson  plot.   The  inside 
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and  outside  heat  transfer  coefficients  (h.  and  h  )  are  then 
solved  for  as  they  were  in  the  plain-end  reduction  except  as 
modified  by  the  different  geometry.   The  equations  used  are 
listed  in  Appendix  D. 

Equation  (36)  was  verified  by  comparing  the  re- 
sults of  the  copper-nickel  and  aluminum  smooth  tubes  for 
equation  (36)  to  the  results  obtained  using  the  equations 
derived  from  equation  (8). 

c.  Friction  Factor 

The  analysis  for  this  reduction  scheme  is  identi- 
cal to  the  smooth  end  reduction  scheme  except  as  modified  by 
the  different  geometry  used. 

d.  Performance  Criteria 

The  area  ratios  were  not  used  in  the  reduction 
based  on  the  hydraulic  diameter. 

The  Colburn  Analogy  was  used  exactly  as  for  the 
smooth  end  reduction  case  described  earlier,  and  shown  in 
Appendix  D. 

3.   Computer  Program 

A  computer  program  was  developed  to  facilitate  both 
rapid  and  consistent  reduction  of  the  data.   A  flow  diagram 
is  shown  in  Figure  19.   Appendix  H  is  a  listing  of  the  pro- 
gram.  The  program  can  essentially  be  divided  into  four 
regions. 
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Region  1.  In  this  section,  the  transport  properties 
and  flow  characteristics  are  determined.  The  equations  used 
for  the  transport  properties  were  obtained  from  reference  [22] 

Region  2.   In  this  section,  the  heat  transfer 
coefficients  are  found. 

Region  3.   The  enhanced  section  pressure  drop  and 
resulting  friction  factor  are  found  in  this  section. 

Region  4.   In  this  section,  the  ratio  of  heat  trans- 
fer performance  to  friction  factor  performance  is  found. 
The  area-ratio  equations  were  not  programmed  at  this  time, 
but  could  be  included  at  a  future  date. 

The  program  development  was  such  that  the  data  would 
be  evaluated  for  both  the  smooth  end  diameter  geometry  and 
the  hydraulic  diameter  geometry  each  time  the  program  was 
activated.   This  was  done  using  a  simple  interation  technique 
and  a  series  of  logic  statements.   The  copper-nickel  smooth 
tube  and  the  aluminum  smooth  tube  were  used  to  proof  the 
hydraulic  diameter  geometry  section  of  the  program  as 
mentioned  above. 
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IV.   RESULTS  AND  DISCUSSION 

A.   INTRODUCTION 

Table  3  lists  the  various  runs  made  and  the  corresponding 
tubes  used  during  these  tests.   Tables  4  through  12  contain 
all  the  data  used  to  evaluate  the  performance  of  the  enhanced 
and  smooth  tubes. 

It  is  seen  in  Table  3  that  several  more  runs  were  made 
for  tube  1  and  tube  4  than  the  other  tubes.   Tube  1  was  the 
first  tube  tested  and,  therefore,  the  procedures  for  testing 
the  enhanced  tubes  were  evaluated  using  this  tube.   Run  3  on 
11  October  1977  was  conducted  with  a  twisted  tape  installed 
at  the  tube  inlet.   These  tapes  were  provided  with  each  tube 
and  were  manufactured  such  that  they  had  the  same  helix  angle 
as  their  respective  tube.   The  results  of  this  run  can  be 
found  in  Table  13.   Basically,  the  pressure  drop  increased  by 
about  15  percent  without  a  corresponding  increase  in  the  heat 
transfer  performance  as  compared  to  the  results  of  run  2.   It 
was  originally  intended  to  include  the  twisted  tape  tests 
within  the  scope  of  this  experiment,  but  due  to  time  constraints 
further  testing  using  these  tapes  was  not  done. 

Run  7  on  7  November  1977  was  conducted  solely  for  the  pur- 
post  of  making  a  movie  of  the  alternate  flooding  and  draining 
between  the  flutes  of  the  45°  HA  enhanced  tube.   The  movie  was 
filmed  such  that  views  of  the  alternate  flooding  and  draining 
would  be  seen  as  viewed  from  the  top  as  well  as  from  the 
bottom  of  the  tube. 
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Figure  20  is  a  sequence  of  six  frames  from  this  movie. 
The  arrow  on  the  top  left  frame  points  to  a  short  length  of 
reflected  light.   In  the  next  frame  down,  the  arrow  points 
to  the  same  spot,  but  the  reflected  light  now  appears  to  be 
about  twice  as  long  as  it  was  in  the  first  frame.   In  the 
third  frame  down,  it  is  seen  that  the  length  of  light  is 
about  the  same  as  it  was  in  the  top  frame.   The  sequence  of 
three  frames  to  the  right  show  the  same  type  of  behavior, 
but  in  a  different  spot  on  the  periphery  of  the  tube.   This 
figure  demonstrates  how  the  condensate  builds  up  in  the 
troughs  to  a  certain  level  and  then  drains  off  as  drops  fall 
off  the  bottom  of  the  tube.   The  set  of  three  frames  on  the 
right  demonstrates  how  this  flooding  and  draining  shifts  from 
trough  to  trough.   As  the  movie  is  viewed,  the  periodic  nature 
of  this  alternate  flooding  and  draining  is  evident/ 

Run  11  was  the  repeatability  run  for  the  45°  HA  tube. 
The  results  of  this  run  were  compared  to  the  results  of  run 
2.   As  can  be  seen  in  Tables  14  and  15,  the  results  are  well 
within  the  experimental  error  band.   The  Sieder  Tate  constant 
for  run  2  is  about  0.081  while  in  run  11  it  is  about  0.078, 
a  difference  of  less  than  4  percent. 

Both  tube  2  and  tube  3  results  were  repeatable  as  can  be 
seen  in  Tables  16  and  17  for  the  30°  HA  tube  and  Tables  18  and 
19  for  the  60°  HA  tube.   Run  10  on  5  January  1978  was  a  poor 
run;  therefore,  the  results  of  this  run  are  not  included. 

Run  4,  taken  with  the  copper-nickel  smooth  tube,  was  con- 
ducted to  insure  that  this  experimenter  could  duplicate  other 
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experimenters'  results  for  a  smooth  tube.   As  can  be  seen  in 
Table  20,  the  Sieder  Tate  constant  was  about  .025  as  compared 
to  .027,  the  value  most  often  found  in  the  literature. 

Runs  5  and  6  were  conducted  to  determine  the  effect  of 
the  wall  thermocouple  location  on  the  wall  temperature  record- 
ed.  It  was  determined  that  the  circumferential  location  of 
the  thermocouple  had  a  profound  impact  on  the  value  of  the 
wall  temperature.   For  this  reason,  all  of  the  data  runs  were 
conducted  with  the  wall  thermocouple  located  on  the  bottom  of 
the  tube,  thus  insuring  that  the  location  of  the  wall  thermo- 
couple would  not  introduce  unaccountable  differences  in  the 
resul ts  . 

Since  a  linear  regression  subroutine  was  used  to  obtain 
the  slope  for  the  Wilson  plot,  the  heat  transfer  information 
obtained  was  very   much  dependent  on  how  well  the  linear  re- 
gression program  could  fit  the  data.   Figure  21  shows  that 
the  generated  curve  fit  the  data  \/ery   well. 

B.   RESULTS  BASED  ON  FLOW  RATES  USING  1.48  GPM  AND  18.8  GPM 
ROTAMETERS 

The  1.48  gpm  rotameter  was  installed  so  that  more  points 
at  low  cooling  water  velocities  could  be  included  in  the 
analysis.   It  was  installed  after  the  initial  data  runs  were 
completed,  since  it  was  felt  that  information  in  the  range 
of  velocities  from  1  to  4  ft/sec  would  be  useful  when  evaluat' 
ing  the  performance  of  the  augmented  tubes.   This  required 
that  the  results  for  velocities  over  the  range  of  both 
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rotameters  be  compared  to  the  results  for  velocities  over  the 
range  of  just  the  18.8  gpm  rotameter.   The  heat  transfer  re- 
sults for  the  range  of  flow  velocities  for  both  rotameters 
did  not  compare  well  with  the  heat  transfer  results  obtained 
for  the  flow  velocities  for  just  the  18.8  gpm  rotameter.   It 
was  found  however  that  by  basing  the  results  of  the  later  runs 
on  the  flow  velocities  over  the  range  of  the  18.8  gpm  rota- 
meter only,  repeatability  of  the  results  for  the  initial  data 
run  was  possible. 

After  a  careful  examination  of  the  results,  it  was  deter- 
mined that  the  reason  for  this  anomaly  was  due  to  variations 
in  the  outside  heat  transfer  coefficient  at  low  velocities. 
Since  at  low  flow  velocities  the  average  cooling  water  tem- 
perature increases,  the  tube  wall  temperature  also  increases 
which  in  turn  causes  the  outside  heat  transfer  coefficient 
to  vary  significantly  as  can  be  seen  from  Nusselt's 
equation  [15]: 
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The  overall  heat  transfer  coefficient,  based  on  the  outside 
area,  is  a  function  of  h . ,  h  ,  and  Rw  as  discussed  earlier 

and  as  seen  in  equation  (7)    which  is  re-written  below: 
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As  was  shown  in  Chapter  III,  equation  (7)  is  reduced  to 
the  form  of  a  straight  line  where  h.  is  a  function  of  the 
Reynolds  number  and  thus  the  water  velocity,  assuming  both 
the  wall  resistance  and  h   remain  constant  over  the  range  of 
water  velocities  used.   The  fact  that  h  ,  as  shown  by  equa- 
tion (37)  above,  may  not  be  constant  at  low  water  velocities 
violates  one  of  the  basic  requirements  of  the  Wilson  Plot 
technique.   Since,  at  low  velocities,  it  is  impossible  to 
keep  the  variation  in  h   small,  the  Wilson  Plot  for  the  flow 
velocities  of  the  1.48  gpm  and  18.8  gpm  rotameters  together 
will  yield  a  significantly  different  value  for  the  Sieder 
Tate  constant  than  for  the  velocities  over  one  rotameter's 
range.   This  is  seen  in  Figure  22.   By  superimposing  a 
straight  line  over  the  1.48  gpm  rotameter  points,  it  is  ob- 
vious that  the  slope  of  the  straight  line  would  be  greater 
than  that  of  the  straight  line  obtained  from  the  linear  re- 
gression subroutine  [18].   On  the  other  hand,  if  a  straight 
line  is  superimposed  over  the  18.8  gpm  rotameter  points,  a 
slope  somewhat  less  would  be  realized.   The  smaller  slope 
would  in  turn  result  in  a  higher  Sieder  Tate  constant  as  seen 
in  equation  C 1  0  c  )  .   This  is  precisely  what  happens  when  the 
1.48  gpm  rotameter  points  are  not  used  in  evaluating  the  data 

Since  the  initial  data  runs  were  taken  without  the  1.48 
gpm  rotameter  points,  and  since  the  pressure  drop  information 
is  felt  to  be  more  reliable  at  the  higher  velocities  through 
the  test  tube,  the  results  that  follow  therefore,  are  pre- 
sented based  on  the  18.8  gpm  rotameter  points  only  unless 
specified  otherwise. 
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PRESSURE  DROP  RESULTS  USING  THE  PRESSURE  TAP  TUBES 


The  General  Atomic  Co.  was  successful  in  drilling  three 
0.024  inch  diameter  holes  into  the  test  section  of  three 
spirally  fluted  tubes  similar  to  the  heat  transfer  test  tubes. 
(See  Figures  14  and  15.)   These  tubes  were  connected  to  the 
pressure  drop  measurement  system  independently,  and  the  pres- 
sure drops  were  measured  at  different  water  flow  velocities 
through  these  tubes.   The  pressure  drop  measurements  are 
summarized  in  Tables  21,  22,  and  23.   Figure  23  is  a  plot  of 
the  pressure  drop  for  60  percent  flow  through  the  18.8  gpm 
rotameter  for  each  of  the  three  different  General  Atomic  tubes 
Note  the  linear  variation  of  the  pressure  drop  along  the 
length  of  the  tubes.   Using  this  pressure  drop  information, 
it  was  possible  to  compute  the  friction  factor  inside  the  test 
section  without  having  to  apply  any  corrections.   It  was  also 
possible  to  examine  the  actual  expansion  and  contraction  los- 
ses caused  by  the  abrupt  geometry  change  in  the  test  tube. 

The  expansion  and  contraction  losses  were  obtained  by 
computing  the  appropriate  pressure  drops  using  equations  (38) 
and  (39)  below: 
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pressure  tap  tubes  as  seen  in  Figure  15,  and  where 
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:   is  the  pressure  drop  across  the  smooth 
end  of  the  tube  between  the  hydraulic 
diameter  section  and  tap  2  or  tap  6. 
Now  the  abrupt  exit  and  entrance  coefficient   can  be  deter- 
mined using  equations  (42)  and  (43): 
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These  values  were  found  at  specific  flow  velocities  and  then 

an  average  of  these  values  was  used  in  reducing  the  overall 

pressure  drop  to  the  pressure  drop  in  the  three  foot  enhanced 

section.   Tables  24,  25,  and  26  summarize  this  procedure. 

Knowing  K  „  and  K  and  using  equation  (12)  the  friction 
3   en      e         3i 

factor  inside  the  enhanced  section  was  determined  as  outlined 
in  Chapter  III.  Tables  27,  28,  and  29  summarize  the  friction 
factors  for  the  PT  tubes.  See  Appendix  H,  sample  calculations 
for  PT  tubes,  for  an  in-depth  demonstration  of  the  procedure. 
As  can  be  seen  when  comparing  the  friction  factors  found 
in  Table  14  to  those  found  in  Table  28,  the  friction  factors 
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found  by  reducing  the  pressure  drop  using  the  values  of  K 


en 


and  K  as  found  in  [20],  are  significantly  different  than  the 

friction  factors  found  by  reducing  the  pressure  drop  using 

the  values  of  K  „  and  K  found  in  this  experiment, 
en      e  r 

The  pressure  drop  across  the  PT  tubes  was  significantly 
higher  than  the  pressure  drop  across  the  corresponding  tubes 
manufactured  for  the  heat  transfer  measurements.   The  most 
obvious  reason  for  this  difference  in  pressure  drop  was  due 
to  the  differences  in  the  internal  geometry  of  the  pressure 
tap  tubes  as  compared  to  the  heat  transfer  tubes.   The  PT 
tubes  were  not  from  the  same  stock  that  the  actual  test  tubes 
were  from,  but  in  fact  they  were  tubes  that  had  been  manu- 
factured earlier  in  order  to  test  out  the  manufacturing 
technique.   Since  there  is  no  convenient  way  of  verifying 
the  geometry  of  the  PT  tubes,  the  results  presented  hereafter 
are  based  upon  the  pressure  drop  measurements  taken  with  the 
heat  transfer  tubes,  using  the  data  reduction  scheme  mentioned 
earlier.   This  procedure  was  used  knowing  that  there  is  an 
error  in  the  expansion  and  contraction  loss  factors. 

D.   RESULTS  BASED  ON  THE  SMOOTH  END  DIAMETER,  D. 

In  as  much  as  the  primary  purpose  of  this  experiment  was 
to  evaluate  the  performance  of  the  augmented-condenser  tubes 
manufactured  by  General  Atomic,  composite  plots  comparing  the 
various  performance  parameters  of  these  tubes  and  smooth  tubes 
are  provided.   Although  it  would  be  highly  desirable  to  com- 
pare these  special  tubes  with  other  augmented  tubes  of 
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different  geometries,  comparison  with  existing  results  was 
avoided  in  most  cases  due  to  differences  in  experimental  con- 
ditions.  The  reason  for  this  is  due  to  the  uniqueness  of 
this  test  facility  as  well  as  the  fact  that  there  are  prob- 
able differences  in  the  data  reduction  schemes. 

1.   Heat  Transfer  Coefficients 

The  first  comparison  is  in  Figure  24,  which  shows 
the  corrected  overall  heat  transfer  coefficient  versus  mass 
flow  rate.   As  expected,  the  enhanced  tubes'  overall  heat 
transfer  coefficients  are  superior  to  the  overall  heat  trans- 
fer coefficients  for  the  smooth  tubes.   At  a  cooling  water 
mass  flow  rate  of  1.3  lbm/sec,  it  is  seen  that  the  smooth 
tubes  (both  aluminum  and  copper-nickel)  have  a  corrected 
overall  heat  transfer  coefficient  of  approximately  1100  BTU/ 
hr-ft   °F  and  that  the  45°  HA  enhanced  tube  has  a  coefficient 
of  approximately  1689  BTU/hr-ft  °F.   This  is  a  53  percent 
increase  in  corrected  overall  heat  transfer  coefficient  for 
the  General  Atomic  tube  over  the  smooth  tubes.   If  the  over- 
all heat  transfer  coefficients  of  the  three  augmented  tubes 
are  compared,  it  is  seen  that  the  60°  HA  tube  (Table  3,  Tube  3) 
outperforms  the  other  two  augmented  tubes.   The  reason  for 
this  is  clearly  seen  in  Tables  14,  16,  and  18.   The  outside 
heat  transfer  coefficient  is  higher  for  the  60°  tube,  thus 
the  outside  thermal  resistance  is  smaller.   As  can  be  seen  in 
equation  (7)  this  makes  U  larger.   The  probable  cause  of  this 
is  that  the  steeper  helix  angle  allows  for  better  condensate 
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drainage  off  this  tube  when  compared  to  the  30°  and  45°  tubes. 
This  fact  will  be  discussed  in  more  detail  in  a  following 
section. 

The  corrected  overall  heat  transfer  coefficient  of  the 
copper-nickel  smooth  tube  and  the  aluminum  smooth  tube  compare 
well  up  to  about  1.3  Ibm/sec.   At  higher  mass  flow  rates, 
there  appears  to  be  some  divergence  in  the  data.   It  is  sus- 
pected that  an  oxide  film  builds  up  on  the  exterior  of  the 
aluminum  tube  more  than  on  the  copper-nickel  tube.   This  in 
turn  caused  the  outside  thermal  resistance  of  the  aluminum 
tube  to  be  higher  than  the  corresponding  resistance  in  the 
copper-nickel  tube.   This  can  be  seen  by  referring  to  Table  20 
and  Table  30  where  the  outside  heat  transfer  coefficient  for 
the  aluminum  tube  is  less  than  for  the  copper-nickel  tube. 
The  inside  heat  transfer  coefficients  are,  however , ^comparable 

2.   Pressure  Drop  and  Friction  Factors 

Figure  25  is  a  comparison  of  the  pressure  drops  versus 
mass  flow  rate  for  the  tubes  tested.   As  can  be  seen,  the 
pressure  drops  of  the  enhanced  tubes  increase  at  a  much  faster 
rate  than  the  pressure  drop  for  the  smooth  tubes.   Neither  the 
30°  HA  tube  nor  the  45°  HA  tube  has  as  dramatic  a  pressure 
drop  as  the  60°  HA  tube.   The  reason  for  this  is  easily  seen 
when  referring  to  Figures  12  and  13.   The  flutes  on  the  30   HA 
tube  have  the  most  gentle  departure  from  the  tube  axis  while 
the  flutes  on  the  60°  HA  tube  have  the  most  profound  departure 
Therefore,  the  flutes  on  the  60°  HA  tube  approach  the  geometry 
of  straight  radial  fins  which  merely  protrude  into  the  inside 


56 


of  a  smooth  tube,  causing  significant  turbulent  mixing  to 
occur.   The  flutes  on  the  30°  HA  tube  are  such  that  they  tend 
to  impart  more  swirl  to  the  fluid,  with  less  turbulent  mixing. 

These  facts  are  further  evidenced  by  referring  to 
Figure  26.   The  friction  factor  for  the  60°  HA  tube  at  a 
Reynolds  number  of  40,000  is  approximately  0.057.   The  friction 
factor  for  the  smooth  tube  at  this  Reynolds  number  is  about 
0.0059.   The  friction  factor  of  the  augmented  tube  is  then 
nearly  ten  times  that  of  the  smooth  tube.   The  45°  HA  tube 
and  the  30   HA  tube  do  not  experience  as  great  an  increase  in 
the  friction  factor,  however,  the  rise  is,  nevertheless, 
significant.   As  expected  the  friction  factor  for  the  30°  HA 
tube  has  the  smallest  value  among  the  enhanced  tubes  for 
Re  >  10,000.   Below  this  Reynolds  number,  the  friction  factor 
for  the  30°  HA  tube  is  greater  than  for  the  45°  HA  tube.   The 
reason  for  this  was  felt  to  be  caused  by  inaccuracies  in  read- 
ing the  pressure  drop  at  low  flows. 

3.   Performance  Criteria 


In  1883,  Reynolds  [23]  mathematically  expressed,  for 
smooth  tubes,  the  analogy  between  heat  transfer  and  momentum 
transfer  as; 

St   =  \  (44) 

In  1933  Colburn  [24]  extended  this  analogy  to  include  Prandtl 
number  effects: 
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The  Colburn  Analogy  then  provides  a  useful  method  to  compare 
the  performance  of  the  enhanced  tube  to  the  smooth  tube. 
Figure  27  is  such  a  comparison.   It  shows  the  tube  perform- 
ance factor,  2j/f  as  a  function  of  Reynolds  number.   As  can 
be  seen,  the  ratio  2 j / f  is  less  than  one  for  all  three  of  the 
enhanced  tubes,  which  indicates  that  the  performance  of  these 
tubes  is  inferior  to  that  of  a  smooth  tube.   But,  as  will  be 
seen  later,  the  tube  performance  factor  is  highly  dependent 
on  which  tube  diameter  is  used  to  evaluate  the  heat  transfer 
and  friction  performance.   The  most  important  fact  to  glean 
from  this  figure,  then,  is  the  relative  performance  of  the 
three  enhanced  tubes.   The  30   HA  tube  clearly  outperforms 
the  other  two.   The  reason  for  this  can  be  readily  seen  in 
Figures  28  and  26.   Figure  28  shows  that  the  inside  heat 
transfer  performance  is  very    nearly  constant  for  the  three 
enhanced  tubes  while  in  Figure  26,  the  frictional  resistance 
of  the  30°  HA  tube  is  seen  to  be  much  less  than  the  other 
two  tubes. 

An  area  ratio  such  as  defined  by    equation  (22)  in 
Chapter  III  can  be  found  from  Figure  28.   In  this  figure,  the 

Sieder  Tate  correlation  is  used  in  place  of  the  Dittus  Boelter 

-1/3 
correlation  as  used  in  equation  (22).   A  value  of  (Nu  Pr 

[u/u  J     )  is  found  at  the  same  Reynolds  number  for  both  the 
w 

smooth  tube  and  the  enhanced  tubes.   The  area  ratio  is  then 
found  by  taking  a  ratio  of  these  dimensi on! ess  groups.   As 
can  be  seen,  this  area  ratio  would  make  the  enhanced  tubes 
appear  to  be  yery   desirable  for  condenser  use. 
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However,  as  explained  in  Chapter  III,  with  this  area 
ratio,  the  additional  frictional  resistance  of  the  enhanced 
tubes  is  ignored.   By  evaluating  the  area  ratio  for  a  con- 
stant pumping  power,  these  additional  frictional  effects  are 
included  in  the  results  as  shown  in  Figure  29. 

The  results  in  this  figure,  however,  are  for  an  exter- 
nal resistance  equal  to  zero.   As  shown  earlier,  area  ratios 
for  a  non-zero  external  resistance  could  also  be  found,  and 
these  results  are  plotted  in  Figure  30. 

When  Figures  29  and  30  are  reviewed,  the  30°  HA  tube 
appears  to  have  the  best  overall  performance.   The  area  ratio 
for  this  tube  is  under  one  over  the  entire  range  of  Reynolds 
numbers  used.   In  the  range  of  Reynolds  numbers  that  are  of 
interest  to  the  condenser  designer  (Re    <    50,000)  the  highest 
area  ratio  is  for  the  45°  HA  tube  and  is  approximately  0.67. 
There  are  three  facts  of  interest  to  note  here.   First,  it 
appears  that  any  of  the  three  tubes  tested  would  outperform 
the  smooth  tube.   Secondly,  the  area  ratio  for  the  60°  HA  tube 
is  less  than  the  area  ratio  for  the  45°  HA  tube  in  Figure  30, 
while  the  opposite  is  true  in  Figure  29.   The  reason  for  this 
is  felt  to  be  caused  by  the  higher  external  resistance  in  the 
45°  HA  tube  than  in  the  60°  HA  tube.   This  is  seen  in  Table  14 
and  Table  18  wherein  h   is  found  to  be  on  the  average, 
200  BTU/hr'ft2  °F  larger  for  the  60°  HA  tube  than  for  the  45° 
HA  tube.   It  is  noted  in  Figure  29  that  the  area  ratios  are 
very   nearly  constant.   This  interesting  fact  apparently  sig- 
nifies that  the  ratio  of  heat  transfer  performance  to  hydro- 
dynamic  performance  is  constant,  much  the  same  as  for  the 
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smooth  tube.   Finally,  when  comparing  Figures  29  and  30,  the 

taking  into  account  of  Rext  has  a  significant  effect  on  the 

results.   The  area  ratio,  as  should  be  expected,  will  increase 

when  the  wall  resistance  is  taken  into  account. 

The  uppermost  curve  and  the  lowest  curve  on  Figures 

29  and  30  are  reproduced  from  Figure  9,  Bergles  [3],  and  are 

for  an  assummed  R   ,  =  0.   These  curves  were  reproduced  herein 

ex  t  r 

to  indicate  the  relative  performance  of  the  General  Atomic 
tubes  against  the  tubes  included  in  reference  [3]. 

Figures  31  and  32  are  Figures  1  and  2  from  reference 
[3]  with  the  30   HA  tube's  results  plotted  thereon.   In 
Figure  31,  it  is  seen  that  the  30°  HA  tube  compares  very 
favorably  with  the  tubes  included  in  reference  [3].   Figure  32 
shows  the  friction  factor  of  the  30°  HA  tube  to  be  slightly 
above  the  average  of  the  existing  data.   Of  the  tubes  appear- 
ing on  Figure  31  and  Figure  32,  tube  9,  a  repeated  rib  tube 
tested  by  Webb  [25]  apparently  comes  the  closest  to  duplicat- 
ing the  General  Atomic  tube's  performance.   However,  it  must 
be  remembered  when  comparing  these  tubes  that  they  are  all 
tested  under  different  conditions  which  could  affect  their 
performance. 

E.   RESULTS  BASED  ON  THE  HYDRAULIC  DIAMETER,  Dh 

The  major  thrust  of  this  section  is  to  permit  the  compar- 
ison of  the  augmented  tubes  based  on  their  respective  hydraulic 
diameters,  in  lieu  of  comparing  them  based  solely  on  their 
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smooth  end  inside  diameter.   Consequently,  it  is  recognized 
that  this  section  will  be  of  interest  primarily  to  heat 
transfer  experimenters. 

As  mentioned  previously,  the  diameter  used  to  compute 
such  parameters  as  the  Reynolds  number  and  Nusselt  number  is 
extremely  critical  to  the  results.   In  reducing  the  data  ob- 
tained in  this  experiment,  the  hydraulic  diameter  was  found 
from  the  geometry  of  the  enhanced  section.   As  expected,  the 
hydraulic  diameter  was  much  smaller  than  the  plain  end  dia- 
meter.  The  effect  of  this  is  immediately  seen  in  Table  31 
through  Table  36,  the  summaries  of  the  enhanced  section  re- 
sults.  For  example,  for  equal  mass  flow  rates  the  Reynolds 
number  will  decrease  and  the  mass  flow  rate  per  unit  area 
will  increase  when  compared  to  the  smooth  end  diameter  results 
This  in  turn  will  cause  the  other  results  to  Ghange  as  well. 

1.   Heat  Transfer  Results 

The  inside  and  outside  heat  transfer  coefficients  are 
both  smaller  in  value  for  the  results  based  on  the  hydraulic 
diameter  in  comparison  to  the  smooth  end  results.   The  major 
reason  for  this  is  due  to  the  fact  that  the  actual  surface 
areas  of  the  enhanced  sections  are  larger  than  the  surface 
area  at  the  smooth  ends.   As  shown  in  Chapter  III,  the  heat 
transfer  rate  can  be  computed  as: 

Q   =   Un  An  LMTD   =   UQ  AQ  LMTD   .  (2) 

For  a  measured  value  of  Q  and  LMTD,  the  UA  product  must 
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remain  constant.   Using  equation  (32),  it  is  easily  seen  that 
if  A-  and  A  both  increase  when  using  the  hydraulic  diameter 
reduction  scheme,  it  follows  that  the  inside  and  outside  heat 
transfer  coefficients  must  decrease.   In  addition,  as  would 
be  expected,  the  Nusselt  number  and  Stanton  number  also  de- 
crease as  seen  in  the  tabular  results. 

2 .   Friction  Factor 

The  friction  factor  found  using  the  hydraulic  diameter 
is  much  less  than  the  corresponding  friction  factor  using  the 
smooth  end  diameter,  as  seen  when  comparing  Figure  33  with 
Figure  26.   The  reason  for  the  smaller  friction  factor  is 
seen  in  the  following  friction  factor  equation: 


Cpb)(.APTS)(144)(2gc) 


TS 


(46) 


4(LTS/D)  G' 


Since  G,  the  mass  rate  of  flow  per  unit  area,  is  inversely 

proportional  to  diameter  squared,  then  the  friction  factor 

5 
is  proportional  to  D  .   Since  D,  is  less  than  D  .  ,  then  the 

friction  factor  will  decrease  accordingly.   As  can  be  seen 
in  Figure  33,  however,  the  friction  factor  for  the  30   HA  tube 
is  less  than  the  corresponding  friction  factor  for  a  smooth 
tube.   Since  this  is  not  reasonable,  the  hydraulic  diameter 
used  herein  may  be  somewhat  questionable.   The  60  HA  friction 
factor  is  again  yery   high  and  apparently  becomes  asymptotic 
at  a  Reynolds  number  of  about  50,000.   Both  the  friction 
factor  curve  for  the  30°  HA  tube  and  the  friction  factor 
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curve  for  the  45   HA  tube  have  a  similar  shape  to  the  curve 
for  the  smooth  tube  as  seen  in  Figure  33. 

3.   Performance  Criteria 

As  seen  in  Figure  34,  the  tube  performance  factor 
2j/f,  when  using  the  hydraulic  diameter,  increases  signifi- 
cantly for  all  three  helix  angle  tubes  when  compared  to  the 
results  based  on  the  smooth  end  diameter. 

F.   ISOTHERMAL  FRICTION  FACTOR  VERSUS  NON- ISOTHERMAL  FRICTION 

FACTOR 

Pence  [13]  compared  the  isothermal  friction  factor  with 
the  non-isothermal  friction  factor  for  the  smooth  tube.   His 
results  showed  a  negligible  difference.   Isothermal  friction 
factors  were  also  found  for  each  of  the  enhanced  tubes  tested 
in  this  experiment.   The  results  are  summarized  in  Tables  37, 
38,  and  39.   Figure  35  compares  the  two  friction  factors  for 
the  45   HA  tube.   As  can  be  seen,  the  isothermal  friction 
factor  is  higher  than  the  non-isothermal  friction  factor  up 
to  a  Reynolds  number  of  about  40,000.   At  this  value  they 
apparently  merge  and  remain  approximately  the  same. 

The  reason  for  the  lower  friction  factor  with  heat  addi- 
tion at  the  lower  flow  velocities  is  thought  to  be  a  result 
of  the  decrease  in  the  measured  pressure  drop  for  these  runs 
as  compared  to  the  isothermal  runs.   The  reason  for  the  reduced 
pressure  drop  must  mean  that  less  energy  is  expended  to  force 
the  water  through  the  enhanced  tube  under  heated  conditions. 
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Consequently,  there  must  be  some  phenomenon  occurring  within 
the  tube  to  aid  the  flow.   It  is  felt  perhaps  that  secondary 
flows  generated  within  the  flutes  are  enhanced  by  the  heat 
addition.   This  fact  was  speculated  by  General  Atomic  and 
Figure  35  supports  this.   At  the  lower  flow  velocities,  the 
heat  input  per  1 bm  of  water  is  greater  than  at  the  higher 
flow  rates  where  the  water  residence  time  in  the  tube  is 
less.   Thus  it  would  be  expected  that  if  the  heat  input  did 
in  fact  enhance  the  secondary  flow  in  the  flutes,  the  differ- 
ence in  the  isothermal  and  non-isothermal  friction  factors 
would  be  greatest  at  the  lower  velocities,  as  verified  in 
Figure  35. 

G.   TUBE  DRAINAGE  AND  HOW  IT  AFFECTS  hQ 

As  is  yery   evident  when  reviewing  the  summary  of  results 
in  the  tables,  the  outside  heat  transfer  coefficient  h 
decreases  for  the  enhanced  tube  when  compared  to  the  results 
for  the  smooth  Cu/Ni  tube. 

As  noted  earlier,  one  very   important  factor  is  the  oxide 
film  that  apparently  forms  on  the  aluminum  surfaces.   Based 
on  the  differences  in  h   for  the  results  of  run  4  and  run  14, 
the  degradation  could  be  as  high  as  10  percent. 

Perhaps,  though,  a  more  important  factor  to  consider  is 
tube  condensate  drainage.   Because  of  the  exterior  flutes, 
and  the  high  surface  tension  of  the  water,  as  condensation 
takes  place,  these  grooves  tend  to  flood  with  condensate, 
which  drops  off  the  tubes  at  regular  intervals.   As  the 
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space  between  adjacent  flutes  fills  with  condensate,  the 
thermal  resistance  on  the  outside  increases.   It  is  seen  that 
the  60°  HA  tube  has  the  highest  outside  heat  transfer  coef- 
ficient among  the  enhanced  tubes.   The  obvious  reason  for 
this  is  that  the  flutes  are  more  conducive  to  condensate 
drainage  than  either  the  45   or  30°  HA  tubes.   This  is  so 
because  the  angle  is  60°  to  the  tube  axis  which  means  that 
the  effect  of  gravity  on  the  water  in  the  flutes  is  greater 
for  this  tube  than  for  the  others,  thus  providing  for  better 
drainage. 
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V.   CONCLUSIONS 

As  a  result  of  the  above-mentioned  tests,  the  following 
conclusions  are  reached. 

1.  The  overall  heat  transfer  coefficient  of  the  General 
Atomic  Co.  tubes  is  superior  to  that  of  a  smooth  tube.   The 
increase  in  the  overall  coefficient,  when  corrected  for  wall 
resistance,  varied  from  about  1.5  to  1.75  times  the  correspond- 
ing smooth  tube  value  at  the  same  mass  flow  rate  of  cooling 
water. 

2.  Significant  enhancement  was  realized  for  the  inside 
heat  transfer  coefficient  in  the  General  Atomic  tubes  when 
compared  to  the  smooth  tube.   Results  based  on  the  smooth 
end  diameter,  show  an  increase  of  about  3.5  times  the  smooth 
tube's  value.   Results  based  on  hydraulic  diameter  show  an 
increase  of  approximately  2. 

3.  The  outside  heat  transfer  coefficient  of  the  augmented 
tubes  was  approximately  10  to  20  percent  less  than  the  cor- 
responding outside  heat  transfer  coefficient  for  the  smooth 
tube. 

4.  Friction  factors  of  the  enhanced  tubes  were  larger  than 
the  corresponding  smooth  tube  values.   At  a  Reynolds  number  of 
50,000,  based  on  D.,  the  friction  factor  for  the  60°  HA  tube 
increased  by  a  factor  of  10  over  the  value  for  the  smooth  tube. 

5.  The  heat  transfer  coefficients  (h.,  hQ,  Uc)  are  highest 

for  the  60°  HA  tube.   The  friction  factor  is  also  the  highest 
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for  this  tube.   The  increase  in  friction  factor  dominates 

the  increase  in  heat  transfer  coefficients  as  exhibited  by  low 

tube  performance  factors,  2j/f. 

6.  The  30°  HA  has  the  highest  tube  performance  factor 

of  the  three  General  Atomic  enhanced  tubes.   The  tube  perform- 
ance factor  is  about  0.67  for  data  reduced  based  on  D.  and  as 
high  as  2.0  for  data  reduced  based  on  the  hydraulic  diameter. 

7.  Condenser  tube  surface  areas  could  be  reduced  by  up 
to  50  percent  using  the  30   HA  tube.   Significant  material 
savings  are  also  possible  using  the  45°  HA  or  60°  HA  tube. 

8.  Using  the  enhanced  tubes,  the  non-isothermal  friction 
factor  is  less  than  the  isothermal  friction  factor  at  the 
lower  water  velocities.   Perhaps  secondary  flows  generated 

in  the  flutes  are  enhanced  by   heat  addition. 

9.  The  range  of  flow  velocities  over  which  the  Wilson 
plot  can  be  accurately  used  is  limited  because  at  low  flow 
velocities,  it  becomes  virtually  impossible  to  maintain  h 
constant.   The  Wilson  plot  assumes  all  variables  are  constant 
except  for  h  .  . 

10.  The  aluminum  tubes  have  a  higher  outside  thermal 
resistance  due  perhaps  to  an  oxide  film.   The  aluminum  smooth 
tube's  overall  heat  transfer  coefficient  compares  favorably 
with  the  copper-nickel  smooth  tube  up  to  about  60  percent  flow 
in  the  system.   At  this  point  the  two  overall  heat  transfer 
coefficients  diverge  and  the  aluminum  tube  is  about  10  percent 
less  than  the  copper-nickel  tube. 
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11.   The  hydraulic  diameter  used  for  the  enhanced  tubes 
may  be  in  error  since  the  30   HA  tube's  friction  factors 
are  less  than  the  corresponding  friction  factor  for  a  smooth 
tube. 
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VI.   RECOMMENDATIONS 

The  results  of  this  experiment  have  been  encouraging 
and  hopefully  the  greatest  contribution  to  heat  transfer  will 
be  to  create  a  curiosity  that  will  cause  others  to  conduct 
experiments  that  will  begin  to  answer  so  many  of  the  unanswer- 
ed questions  and  of  course,  to  provide  further  experimental 
evidence  to  support  the  results  herein.   In  an  effort  to  offer 
some  ideas  for  further  experiments,  the  following  recommenda- 
tions are  made. 

1.  To  fully  realize  the  potential  of  this  condenser  test 
system  as  presently  configured,  as  many  of  the  augmented  tubes 
available  commercially  as  well  as  newly  conceived  tubes  should 
be  tested  and  ranked  with  other  tubes.   This  will  present  to 
the  designer  of  Naval  condensers,  as  well  as  other  designers, 

a  clear  and  non-proprietary  reference  by  which  to  evaluate 
his  possibilities.   Figures  36  and  37  show  tubes  manufactured 
by  Spiral  Tubing  Corp.  and  Wolverine  Division  of  Universal  Oil 
Products  Co.  respectively,  which  are  examples  of  tubes  that 
can  be  tested  and  ranked  in  the  future. 

2.  The  enhanced  condenser  tubes  used  in  this  study  should 
be  tested  in  a  vertical  orientation.   As  was  discussed  earlier, 
in  a  horizontal  orientation,  the  condensate  fills  the  area 
between  the  flutes  which  has  a  detrimental  effect  on  the  out- 
side heat  transfer  coefficient.   It  is  felt  that  a  vertical 
orientation  would  provide  a  better  outside  heat  transfer 
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coefficient  since  superior  tube  drainage  would  result  in  a 
vertical  orientation.   A  preliminary  design  for  modifying  the 
existing  system  to  allow  a  vertical  orientation  of  the  con- 
denser tube  was  made  by  this  experimenter  and  will  allow  mod- 
ifications at  a  minimal  cost.   Included  in  this  design  are 
also  various  baffle  designs  to  permit  different  angles  for 
steam  entry  into  the  condensing  section  of  the  condenser. 
By  directing  the  steam  at  different  angles  across  the  condenS' 
ing  surface,  an  optimum  entry  angle  may  be  found  which  will 
give  the  highest  outside  heat  transfer  coefficient. 

3.  Tests  of  the  most  promising  tubes  should  be  expanded 
to  include  testing  of  horizontal  banks  of  tubes.   This  type 
of  test  is  needed  to  determine  how  the  condensate  drainage 
of  the  upper  tubes  affects  the  lower  tubes'  heat  transfer 
coefficient. 

4.  Further  hydrodynamic  testing  of  augmented  tubes  is 
required.   The  pressure  drops  and  hence  the  friction  factors 
of  augmented  tubes  cannot  be  accurately  predicted  by  existing 
methods;  therefore,  tests  which  permit  the  exact  measurements 
of  pressure  drops  on  tubes  of  known  dimensions  should  be 
conducted. 

5.  Pressure  drop  measurements  of  the  augmented  tubes 
under  various  degrees  of  heat  addition  to  the  cooling  water 
should  be  accomplished.   This  is  needed  to  determine  how  much 
the  heat  addition  affects  the  secondary  flows  induced  in  the 
flutes,  and  the  resulting  pressure  drop.   Also,  tests  which 
permit  flow  visualization  under  both  heat  transfer  conditions 
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and  isothermal  conditions  would  be  desirable  as  this  type 
of  test  could  possibly  shed  light  on  secondary  flow  enhance- 
ment in  these  tubes. 

6.   The  fouling  characteristics  of  any  tube  under  test 
should  be  determined.   As  discussed  within  this  report,  there 
is  some  evidence  that  the  aluminum  tube  forms  an  oxide  film 
that  increases  the  thermal  resistance.   By  manufacturing 
and  testing  identically  dimensioned  tubes  of  different 
material,  this  effect  could  be  determined.   In  addition  to 
the  outer  film,  the  extent  of  the  fouling  problem  on  the 
cooling  water  side  of  the  tube  has  to  be  determined. 
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VII.     FIGURES 
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FLOW  BAFFLES 


*  CONDENSATE 


Figure  5.   Test  Condenser  Schematic,  Side  View 
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igure  11.   Photograph  of  Forming  Process,  Spirally  Fluted  Tubes 


82 


^T-\  ''.'■'•-'        *      ■■        ':y,^\  ' 


83 


,    ■ 


.  ;         ■■ 


■  ■.'.'.v. .  .-.■■.■  ■  -      ■.;.;.;• 


o 
o 


to 
Qi 

3 


■o 


o 

LO 


CD 


2 
CD 


A3 

C 

o 

•r- 
■M 
O 
<l) 
V) 
t/> 

o 

s- 
o 


ro 


tsUm, 


*&%^->Z:<>-.  :■-■<>  •  :'•  -.;.V:'.;  ^     :-':  .'■.:.,.    :"'";:< v;^.*. 


o 

o 

CO 


84 


85 


«J3 

I— 
CL. 


/v  y\ 


in 

i 


H 


X 


LU 


00 


CO 


^ 

C>|c\J 

lo|c\j 

r—  |<«0 

CO 

r-|oo 

< 

r— |0O 

c\j|co 

' —  l« — 

1— 

C\J 

vo 

LT) 

CO 

00 

CL 

1 

1 

CO 

r-^|co 

r-.|oo 

r-|CvJ 

L0|00 

C\J 

co|oo 

1 

in 

m 

CO 

00 

^— 

co 

T— 

r— 

CO 

Q. 

uoko 

1 

H- 

—  |<NJ 

'—]'— 

H- 

LT) 

VO 

CO 

00 

Q. 

•ei- 

CO 

c 

c 

c 

c 

c 

<U 

J3 

3 

h— 

<c 

en 

o 

Q 

LU 

<L> 


c 

«5 


i. 

3 
«/» 
(/> 

<1J 

s- 
o. 


o 
o 

■M 

E 
d) 

-c 
u 

oo 

LO 


<D 
J- 
3 
CD 


X 


csj 


X 


Nj^ 


86 


o 
-o 

CO 


s_ 

3 

-o 

<D 
U 
O 

S- 

c_ 

M- 
O 

c 
o 

•I— 

-!-> 
03 
-t-> 

C 

<v 
to 

<u 
s_ 

Q. 
CD 


E 
o 

CO 


vo 


i- 
en 


87 


s. 

CD 
■M 
<U 

E 

fO 

s- 

Q- 

+J 
T3 


i- 

■o 

•r* 
C/l 


-a 
cu 


cu 
i_ 

T3 
<U 
U 

o 

Q. 

4- 
O 

C 

o 

•r* 
■U 
<0 

4-> 

C 

a> 

CO 
0) 

'_ 
a. 
a> 


■4-> 


<U 


0> 


88 


-o 

OJ 

to 

• 

=) 

( 

-C 

<u 

s. 

T3 

3 

C 

T3 

<o 

<D 

a 

•r 

o 

J= 

S- 

a_ 

•> 

+-> 

<+- 

C 

o 

JO 

+-> 

c 

10 

o 

c 

•r— 

o 

+-> 

<_> 

«a 

4-> 

CD 

C 

+-> 

O! 

»o 

10 

t— 

<u 

s_ 

s- 

a. 

<v 

<u 

-a 

q: 

<u 

•r" 

o 

oo 

•r— 

4-> 

(0 

E 

<u 

.c 

CJ 

C/5 

00 


O) 


89 


u.   k.    P.     cp 


»(U 


»e(I) 


0(1) 


■  »r 


..,<■>.  !iiiii2 


»t,(D 


«,  T     -  Tc. 


T(I) 


X 


r 

c(l) 

«n 

°N 

3  I 

P     H 

ell)  •   t^7 


0.40 

rl 

««,llp 

1 

1 

'\ 

_    S.04GU)2 

r, 

'    D1 

*       «e2(l) 


,-75 


«  0(1)'  L 


'if-Wh^ 


-> 


v^».0.8,r,/^ 


«  C(l)?  I. 


h,.^..°-V»i^4 


•  •  1       >       1 


o"nir — rr — r- 


IT- 


< 

1                                                                                            i 

LNR£G 

Pr""-<(.'    )°" 


J       •       St  Pr 


tpp   -  Ls-1 
f 


Figure  19.   Computer  Program  Flow  Diagram. 
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Figure  31.   Heat  Transfer  Data  For  Tubes  with 
Various  Types  of  Roughness  versus 
Reynolds  Number. 
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Various  Types  of  Roughness  Versus 
Reynolds  Number. 
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TABLES 


Channel 

Number 

Location 

Channel 

Number 

Location 

40 

Tc. 
l 

47 

Tv 

41 

Tc0 

48 

Tv 

42 

Tc0 

49 

Tv 

43 

Tc0 

50 

T 

V 

44 

Tc0 

51 

Tw 

45 

T 

V 

52 

Hotwell 

46 

T 

V 

Table  1  .  Location  of  Stainless  Steel  Sheathed 
Copper  Constantan  Thermocouples. 
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Oannel 
umber 


Location 


Channel 
Number 


Location 


Hot  Well 
Feedwater  Tank 
Condenser  Window 

Tc„ 


Condensate  Header 
Tc  into  Cooling  Tower 
Tc  out  of  Cooling  Tower 
Cooling  Tower  Ambient 


Table  2.   Location  of  Teflon  Coated  Copper  Constantan  Thermocouples 
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of 

JO 

FLOW 

V°f) 

TW(°R) 

TCi(°F) 

Tc0(°F) 

AP(psi) 

11 

151.4 

574.4 

69.7 

101.8 

.77 

21 

151.2 

565.4 

69.8 

91.1 

2.68 

31 

151.3 

561.0 

69.9 

86.2 

5.50 

41 

151.2 

557.8 

70.3 

83.4 

9.32 

51 

151.3 

555.8 

70.8 

81.8 

13.92 

61 

151.0 

554.2 

71.2 

80.7 

19.56 

71 

150.8 

552.8 

71.6 

79.9 

26.44 

74 

150.8 

552.6 

71.7 

79.8 

29.21 

Table  4.  Raw  Data  for  45°  Tube,  Run  2,  11  OCT  77 


% 

FLOW 

y°F) 

y°R) 

Tc.(°F) 

TC0(°F) 

AP(psi) 

10 

154.7 

573.3 

69.0 

101.30 

1.00 

20 

153.1 

564.3 

68.5 

90.37 

3.14 

30 

151.9 

559.2 

68.3 

84.74 

6.42 

40 

149.7 

556.2 

68.3 

81.50 

11.10 

50 

149.7 

554.2 

68.7 

79.88 

16.61 

60 

150.0 

552.4 

69.2 

78.89 

23.53 

70 

150.3 

551.0 

69.4 

78.04 

31.45 

73 

150.7 

551.1 

69.9 

78.35 

33.68 

• 

Table  5.  Raw  Data  for  45  Tube,  Run  3,  11  OCT  77 

11  2 


of 

FLOW 

T  (°F) 

V 

T   (°R) 

w 

TCi((V) 

Tc/Y) 

AP(psi) 

10 

152.3 

574.7 

66.3 

84.5 

.18 

20 

150.7 

559.4 

66.0 

79.2 

.54 

30 

151.1 

554.4 

65.6 

76.5 

1.11 

40 

151.1 

552.2 

65.8 

74.8 

1.88 

50 

150.7 

551.0 

66.0 

73.8 

2.75 

60 

150.6 

549.0 

66.3 

73.2 

3.82 

70 

149.3 

547.0 

66.5 

72.5 

5.07 

80 

149.1 

546.3 

66.7 

72.2 

6.39 

90 

149.5 

546.1 

66.7 

71.7 

7.91 

100 

149.4 

545.2 

66.9 

71.6 

9.58 

Table  6.  Raw  Data  for  Smooth  Tube,  Run  4,  18  OCT  77 


% 
FLOW 

VCF) 

T  (CR) 

Tc.(°F) 

Tc0(°F) 

AP(psi) 

10 

149.9 

581.6 

64.2 

99.4 

.81 

15 

150.4 

576.3 

64.0 

92.1 

1.63 

20 

150.1 

574.2 

64.2 

87.7 

2.43 

25 

150.5 

571.8 

65.1 

85.0 

3.48 

30 

150.9 

570.8 

66.3 

83.7 

4.77 

40 

149.7 

568.8 

67.4 

81.3 

7.98 

50 

148.6 

567.3 

67.6 

79.2 

12.03 

60 

148.7 

566.6 

67.2 

77.4 

16.70 

70  ' 

148.3 

564.6 

67.2 

76.1 

22.23 

30 

148.6 

563.7 

67.2 

75.4 

27.96 

i 

Table  7.  Raw  Data  for  30°  Tube,  Run  8,  28  NOV  77 
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% 

FLOW 

V°f) 

y°R) 

Tc.(°F) 

Tc0(°F) 

AP(psi) 

10 

154.1 

579.6 

70.5 

106.4 

.91 

15 

151.1 

590.7 

69.6 

96.4 

2.13 

20 

149.3 

584.8 

68.9 

90.8 

3.73 

25 

149.0 

580.2 

68.7 

87.8 

5.59 

30 

148.7 

579.8 

68.5 

85.7 

7.75 

40 

148.2 

577.3 

68.5 

81.9 

13.74 

50 

147.9 

576.0 

68.4 

80.3 

21.75 

60 

147.5 

578.2 

68.7 

79.0 

31.97 

65    • 

147.4 

576.4 

69.1 

78.7 

38.27 

Table  8.  Raw  Data  for  60°  Tube,  Run  9,  29  NOV  77 


FLOW 

V°f) 

y°R) 

TCi(°F) 

Tc0(°F) 

AP(psi) 

15 

152.2 

567.02 

69.9 

95.9 

1.75 

20 

151.3 

563.0 

67.2 

90.0 

2.57 

25 

151.4 

560.1 

67.6 

86.9 

3.86 

30 

151.3 

557.8 

68.0 

84.7 

5.32 

40 

151.5 

554.9 

68.1 

81.7 

9.03 

50 

151.5 

552.8 

68.7 

80.1 

13.35 

60 

151.2 

551.5 

69.0 

79.0 

18.95 

70 

151.1 

549.5 

69.4 

78.6 

25.71 

76    ' 

151.2 

549.0 

69.8 

77.9 

30.03 

Table  9.  Raw  Data  for  45°  Tube,  Run  11,  9  JAN  78 
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% 

FLOW 

V°f) 

y°R) 

Tc.(°F) 

Tc0(°F) 

Ap(psi) 

10.5 

149.7 

592.7 

64.7 

101.4 

1.11 

15 

149.8 

592.0 

64.5 

94.6 

2.18 

20 

150.1 

588.9 

64.4 

89.8 

3.82 

25 

149.9 

589.5 

64.7 

86.3 

5.46 

30 

149.9 

588.0 

64.9 

84.0 

7.80 

40 

150.0 

585.5 

64.9 

80.5 

13.56 

50 

149.8 

585.3 

65.3 

78.4 

21.00 

60 

149.8 

584.8 

65.6 

77.0 

30.97 

66     • 

149.7 

583.2 

66.3 

76.7 

38.13 

Table  10.  Raw  Data  for  60°  Tube,  Run  12,  10  JAN  78 


FLOW 

V°f) 

TW(°R) 

Tc-(°F) 

Tc0(°F) 

AP(psi) 

10 

151.6 

585.7 

64.5 

84.3 

.15 

15 

151.0 

580.7 

64.4 

81.4 

.29 

20 

150.5 

575.4 

64.4 

79.3 

.52 

25 

150.2 

572.0 

64.5 

77.3 

.77 

30 

149.9 

566.4 

68.7 

76.8 

1.10 

40 

149.9 

561.8 

64.9 

75.  ^ 

1.84 

50 

149.7 

559.6 

65.3 

74.0 

2.82 

61 

149.7 

558.9 

65.6 

73.4 

3.98 

70     • 

149.5 

557.1 

66.2 

73.1 

5.12 

81 

149.5 

555.6 

66.3 

72.5 

6.69 

89 

149.6 

555.1 

66.7 

72.5 

7.88 

Table  11.  Raw  Data  for  Smooth  Al  Tube,  Run  14,  14  JAN  78 
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% 

FLOW 

V°f) 

TW(°R) 

Tc.(°F) 

TC0(°F) 

AP(psi) 

10.5 

151.3 

585.0 

61.8 

97.2 

.72 

15 

152.0 

582.3 

65.1 

92.5 

1.36 

20 

151.2 

579.0 

62.1 

86.2 

2.27 

25 

151.5 

576.9 

64.9 

84.8 

3.36 

30 

150.8 

575.4 

62.1 

80.4 

4.71 

40.25 

151.2 

573.1 

62.4 

77.1 

7.82 

51 

151.4 

571.7 

63.3 

75.5 

11.90 

61 

151.3 

570.4 

64.2 

74.6 

16.49 

70     • 

151.1 

570.4 

64.4 

73.9 

21.45 

80 

151.1 

568.3 

64.9 

73.4 

27.55 

Table  12.  Raw  Data  for  30°  Tube,  Run  15,  24  JAN  78 
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%  Flow 

10 

30 

50 

70 

Taps 

C'Hg) 

("Hg  ) 

("Hg) 

("Hg) 

1-7 

2.0 

13.2 

34.2 

64.5 

2-6 

1.9 

13.1 

34.0 

64.4 

2-4 

.7 

4.6 

11.6 

21.4 

2-3 

0 

.3 

.8 

1.4 

3-5 

1.4 

9.2 

23.6 

44.5 

3-4 

.7 

4.6 

11.0 

20.2 

4-5 

.7 

4.9 

12.7 

24.5 

5-6 

.5 

3.7 

10.0 

19.0 

6-4 

1.2 

8.4 

22.6 

43.2 

Table  21.     Summary  of  Pressure  Drops     PT  -  1 ,  30 
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%  Flow 

10 

20 

30 

40 

50 

60 

TAPS 

("Hg) 

("   Hg) 

("  Hg) 

("Hg) 

("Hg) 

("Hg) 

1-7 

2.3 

7.7 

17 

28.8 

43.6 

61.5 

2-6 

2.2 

7.5 

16.8 

28.3 

43.03 

60.5 

2-4 

.9 

3.1 

6.4 

11.5 

17.6 

24.8 

2-3 

0 

.2 

.4 

.8 

1.2 

1.7 

3-5 

1.7 

5.6 

12.5 

21.3 

32.2 

45.3 

3-4 

.9 

2.9 

6.4 

10.9 

16.6 

23.3 

4-5 

.9 

2.8 

6.2 

10.5 

15.7 

22.2 

5-6 

.5 

1.7 

3.8 

6.5 

9.9 

13.9 

6-4 

1.4 

4.5 

9.9 

16.9 

25.3 

35.6 

Table  22.  Summary  of  Pressure  Drops  PT-2,  45 
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%  Flow 

10 

15 

20 

25 

30 

40 

50 

60 

Taps 

("Hg) 

("Hg) 

("Hg) 

("Hg) 

("Hg) 

("Hg) 

("Hg) 

("Hg) 

1-7 

2.4 

5.4 

8.9 

13.9 

19.7 

34.5 

53.5 

78.2 

2-6 

2.3 

5.1 

8.6 

13.6 

19.4 

34.0 

53.0 

77.0 

2-4 

.9 

2.0 

3.5 

5.5 

7.6 

13.9 

21.9 

32.0 

2-3 

0 

0 

0 

.1 

0 

.2 

.4 

.6 

3-5 

1.9 

4.0 

6.7 

10.5 

14.9 

26.2 

41.0 

60.1 

3-4 

.9 

2.1 

3.5 

5.5 

7.7 

13.8 

21.7 

31.7 

4-5 

.8 

2.0 

3.4 

5.2 

7.2 

12.5 

19.5 

28.7 

5-6 

.5 

1.4 

2.2 

3.2 

4.5 

7.8 

11.8 

16.8 

6-4 

1.4 

3.1 

5.3 

5.5 

11.6 

20.1 

31.1 

45.2 

Table  2  3.     Summary  of  Pressure  Drops     PT  -  3,  60 
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% 

Fl  ov\ 


W 


pv 


2gc 


(psi) 


APAc 
(in   Hg) 


APT 
(in  Hg) 


AP         +u 
smooth 

(in  Hg) 


AP 
comp 

(in  Hg) 


en 


AP 
exp 

(in  Hg) 


20 


0.608 


0.555 


0.741 


0.375 


0.0633 


0.534 


0.45 


0.4135 


0.35 


30 


0.608 


1.26 


1.68 


0.837 


0.1345 


1.208 


0.46 


1.108 


0.42 


40 


0.608 


2.18 


2.91 


1.43 


0.2211 


1.899 


0.40 


1.718 


0.37 


50 


0.608 


3.41 


4.55 


2.15 


0.3270 


2.849 


0.39 


2.949     0.40 


60 


0.608 


4.86 


6.49 


3.03 


0.4570 


3.895 


0.37 


4.155 


0.39 


ZKCj 


EKe 


=    0.41 


=    0.39 


Table  24.     Summary  of  K      +  K    determination  PT-2,  45 
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'     /  1  bin  x 
in    (— -) 

0.260 

0.781 

1.30 

1.82 

Res 

13,492 

37,416 

61,178 

83,894 

fs 

0.0073 

0.0060 

0.0054 

0.0051 

r      , 1 bm       \ 
Gs    fePft2) 

610,609 

1,834,142 

3,057,795 

4,282,229 

AP.   (psi) 

0.0082 

0.0607 

0.152 

0.281 

VTS   (0) 

17,100 

51,300 

85,500 

119,700 

VT<~ 

P    2g|    (Psi) 

0.151 

1.36 

3.77 

7.40 

<Kc  +  Vl 

0.40 

0.40 

0.40 

0.40 

<Kc  +  Ke>2 

0.77 

0.77 

0.77 

0.77 

APexp/con/Psi) 

0.0604 

0.544 

1.51 

2.96 

AP       ,        (psi) 
exp/coru  K 

0.116 

1.0472 

2.903 

5.70 

AP2_6(psi) 

0.84 

5.96 

15.47 

29.31 

APTs^(psi) 

0.7714 

5.355 

13.808 

26.069 

AP^^psi) 

0.7758 

4.852 

12.415 

23.329 

P       / 1 bm       \ 
hr-ft 

1,060,726 

3,186,201 

5,311,884 

7,438,923 

fTS 

0.0092 

0.00707 

0.00656 

0.00631 

fTS2 

0.00852 

0.00641 

0.00589 

0.00565 

AP3_5(psi) 

0.637 

4.18 

10.72 

20.25 

f3-5 

0.00867 

0.00632 

0.00582 

0.00561 

Table  27.  Friction  Factor  Summary,  PT-l,  30 
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ni    (sec") 

.530 

.800 

1.050 

1.310 

1.570 

Res 

27,419 

40,176 

52,086 

64,696 

75,074 

fs 

.0063 

.0059 

.0056 

.0053 

.0052 

n       , 1 bm       x 
Gs    4ir'ft2) 

1,274,654 

1,920,332 

!  2,526J311  3,157,827 

3,770,897 
.2081 

AP$  (psi) 

.0288 

.0612 

.1006 

.1488 

VTS  <n7> 

32,796 

49,428 

64,980 

81,216 

96,984 

P  2gc   (psi) 

.555 

1.26 

2.18 

3.41 

4.86 

<Kc  +  Ke>] 

.33 

.33 

.33 

.33 

.33 

<Kc  +  Ke>2 

.80 

.80 

.80 

.80 

.80 

APexp/con/Psi) 

.183 

.416 

.719 

1.13 

1.60 

APexp/con2(Psi) 

.444 

1.01 

1.74 

2.73 

3.89 

AP2_6(psi) 

3.39 

7.65 

12.88 

19.58 

27.53 

APTS  (psi) 

3.19 

7.17 

12.06 

18.30 

25.72 

APTS   (psi) 

2.92 

6.58 

11.04 

16.70 

23.43 

GTS   (            ?' 
15     hr-fr 

2,036,156 

3  £67,574 

4,035,097 

5,044,372 

6,023,701 

'«, 

0.0115 

0.0114 

0.111 

0.0108 

0.0106 

fTS9 

0.0106 

0.0105 

0.0101 

0.0098 

0.0097 

AP3_5(psi) 

2.55 

5.69 

9.69 

14.63 

20.62 

f3-5 

0.0108 

0.0107 

0.0103 

0.0099 

0.0097 

Table  28-  Friction  Factor  Summary,  PT-2,  45 
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,  1  tun  v 

'•'    lsec  ' 

*• 

f 

b". 

r      ,  1  bin       * 

Ai's   (psl) 

vrs  <hF> 

V2 

P-fgC-(PSl) 

.260 


.390 


.520 


.650 


.780 


1.04 


1.30 


1.56 


14^590  20,516  26,384   32,366   38,331  50,009    61,313   73,799 
.00719  .0066 


.0062 


,0061 


.0059    .0056 


.0054 


.0052 


(Kc  *  Ke)2 


'''''cxp/con^01) 
hi'      ,        (psi) 


609,807  915,703  1,221 ,674  1,527,559  1,833 , 453  2,445,4343,057,3183,669 , 066 
.0073 


.0151 


.0252    .0388 


.0540    .0912    .1375 


1907 


14,976     22,464  29,952         37,440     44,928       59,904       74,880       89,856 

4.17 


.261 


,463 


,724 


1.04 


1.85 


0.26 


0.26 


0.26 


0.26 


0.26 


2.89 
0.26 


1.03         1.03 


1.03 


1.03 


1.03 


1.03 


.030 
.119 


.068 


.269 


,120 


188 


,270 


.481 


477 


,746 


1.071 


1.906 


2-6 


(psi 


(psi) 


TS 


&PTSi.(psi) 


G       (- 
°TS   [ 


I  bin 


hr-ft 


2' 


TS- 


TS, 


AiJ3.5(psi) 


2.3 


3.91 


6.19 


8.83 


15.47 


2.217       3.765         5.963         8.506         14.90 


2.016       3.408         5.405         7.705         13.47 


1.03 

.751 
2.977 

24.12 
23.23 
21.01 


0.26 

1.03 

1.084 

4.295 

35.04 
33.77 
30.55 


927,828  1^93,2531,858,792  2,324,200  2,789,622  3,720,759  4,651,749    5,582,531 

.0196 

.0177 

1.92         18.66         27.35 


.0207       .0197 


.0188       .0179 


1.82 


.0200 


.0181 


4.78 


.0198 
.0179 


.0195 
.0176 


.0194 
.0176 


6.76 


f. 


T-5 


0191        .0179 


.0180 


,0177 


.0175 


0175 


0178 


Table  29.  friction  Factor  Summary,  PT-3,60 
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APPENDIX  A 
CALIBRATION  PROCEDURES 

1  .   Rotameters 

The  three  rotameters  installed  in  the  test  condenser  and 
by-pass  line  were  calibrated  by  the  use  of  a  scale  and  a  stop 
watch.   The  procedure  that  was  used  is  as  follows: 

a.  Establish  flow  through  the  rotameter  at  a  specified 
percent  1 evel  . 

b.  Record  weight  of  container  at  time  zero. 

c.  Record  weight  of  container  at  time  zero. 

d.  Determine  flow  rate  for  specified  level  by 
dividing  the  weight  gain  by  the  elapsed  time. 

2.   Thermocouples 

All  thermocouples  used  in  the  test  apparatus   were  cali- 
brated in  a  silicone  oil  bath  against  a  platinum  resistance 
thermometer.   The  output  of  the  platinum  resistance  ther- 
mometer was  determined  by  the  use  of  a  Wheatstone  bridge  and 
is  accurate  to  the  nearest  0.02°F  (0.01°C).   The  Autodata 
Nine  Recorder,  the  output  device  for  all  sheathed  thermo- 
couples, is  accurate  to  the  nearest  0.2  F  (0.1  C)  and  all 
sheathed  thermocouples  were  found  to  record  the  same  temper- 
ature that  was  indicated  by  the  platinum  standard. 
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3.   Pressure  Transducer 

Calibration  of  the  absolute  pressure  transducer  was  done 
by  using  a  mercury  manometer.   The  amplifier  for  the  trans- 
ducer had  an  excitation  voltage  of  16.2  volts.   A  second 
voltage  setting  was  used  to  set  the  scale  of  the  output  so 
that  the  output  corresponded  to  centimeters  of  mercury.   The 
calibration  showed  a  linear  relationship  between  the  trans- 
ducer and  the  actual  pressure.   The  output  of  the  transducer 
when  multiplied  by  ten  is  equal  to  the  pressure  in  centimeters 
of  mercury. 
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APPENDIX  B 
PROCEDURES  FOR  PREPARING  AND  INSTALLING  THE  ALUMINUM  TUBES 

The  enhanced  section  of  the  tubes  manufactured  by  General 
Atomic  Company  had  a  larger  outside  diameter  than  the  corres- 
ponding smooth  end.   This  fact  necessitated  that  a  special 
procedure  be  developed  to  allow  installation  of  these  tubes 
with  the  wall  thermocouple  installed.   The  smooth  aluminum 
tube  was  prepared  in  a  manner  similar  to  the  aluminum  augmented 
tubes,  except  for  the  special  installation  procedures  mentioned 
above.   Therefore,  these  procedures  also  apply  to  the  smooth 
aluminum  tube, except  step  A.  1  is  ignored  and  all  of  step  C  is 
disregarded. 

A.   PREPARATION  PROCEDURES 

1.  File  the  plastic  thermocouple  connector  to  a  size  that 
will  allow  the  connector  to  pass  through  the  swagelock  fitting. 

2.  Cut  a  4-inch  long  groove  approximately  0.025  inch  deep 
in  the  wall  of  the  tube.   The  groove  is  cut  from  the  edge  of 
the  enhanced  section  back  into  the  plain  end  section  at  one 
end . 

3.  Install  the  thermocouple  bead  at  approximately  the 
middle  of  the  enhanced  section.   Insure  that  the  bead  location 
is  orientated   90°    to  the  0.025  inch  deep  groove  and  that 
the  bead  is  in  contact  with  the  tube.   This  insures  that  the 
thermocouple  bead  will  be  at  the  same  location  on  all  tubes. 
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Use  an  appropriate  high  temperature  cement,  such  as  DEVCON, 
to  install  the  bead. 

4.  Cement  the  thermocouple  bead  into  the  0.025  inch 
deep  groove. 

5.  After  the  cement  dries  ,  insure  that  the  excess  ce- 
ment is  removed  by  scraping  with  a  pocket  knife  and  filing 
smooth  so  as  to  allow  a  smooth  fit  in  the  swagelock  fitting. 

B.   TUBE  CLEANING  PROCEDURE 

Prior  to  each  run,  the  aluminum  test  tube  must  be  cleaned 
thoroughly.   The  procedure  used  is  as  follows: 

1.  Lightly  brush  the  exterior  of  the  tubes  with  a 
wire  brush. 

2.  Fill  the  flood  tube  2/3  full  (The  flood  tube  is  a 
1-inch  ID,  48  inch  long  plexiglass  tube  with  a  1-inch  thick 
piece  of  flat  plexiglass  glued  to  one  end.)   with  a  10  percent 
solution  of  sodium  hydroxide  (NaOH). 

3.  Insert  the  test  tube  into  the  flood  tube. 

4.  Using  a  3/4  inch  test  tube  brush,  brush  the  inside 
surface  of  the  test  tube.   The  brush  must  be  inserted  with  a 
twisting  motion  to  insure  that  the  flutes  are  being  cleaned. 

5.  Leave  the  test  tube  in  the  flood  tube  until  a  foaming 
action  is  observed  over  the  entire  outside  surface  and  the 
solution  has  bubbled  up  through  the  inside  of  the  test  tube. 
This  should  take  about  two  minutes. 
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6.   After  removing  the  test  tube  from  the  flood  tube, 
thoroughly  rinse  the  inside  and  the  outside  with  tap  water. 
The  tube  should  have  a  dull  silvery  appearance  and  should  be 
well  wet  by  the  water. 

C.   INSTALLATION 


The  fluted  aluminum  test  tubes  were  installed  by  following 
the  procedure  listed  below: 

1.  Insertthe  tube  into  the  test  condenser  through  the 
tube  sheet  which  the  thermocouple  lead  will  pass  through. 

2.  Install  the  test  tube  swagelock  fitting  and  seals 
over  the  thermocouple  case  and  the  thermocouple  lead.   Do  not 
allow  the  lead  to  be  seized  by  the  swagelock  fitting  as  this 
will  destroy  the  thermocouple. 

3.  Install  the  upper  dummy  tube's  swagelock  fitting. 

4.  Connect  the  cooling  water  lines  to  the  test  tube. 

5.  Reinsulate  the  test  condenser  and  the  cooling  water 
lines. 
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APPENDIX  C 
OPERATING  PROCEDURES 

1.   Light-off  Procedure 

a.   Boiler  Operation 

(1)  Energize  main  circuit  breaker  located  in  power 
panel  P-2. 

(2)  Turn  key  switch  on--located  on  right  side  of 
main  control  board. 

(3)  Energize  circuit  breaker  on  left  side  of  main 
control  panel  by  depressing  start. 

(4)  Energize  individual  circuit  breakers  on  left 
side  of  main  control  panel.   The  following  list 
identifies  each  circuit  breaker: 

(a)  Feed  pump 

(b)  Outlets 

(c)  Hot  water  heater  (feedwater  tank) 

(d)  Condensate  pump 

( e )  Boiler 

(f)  Cooling  tower 

(g)  Cooling  water  pump  (only  when  using  closed 
cooling  water  system) 

(5)  Insure  water  level  is  up  in  the  feedwater  tank. 
Turn  switch  on  to  energize  heater. 

(6)  Energize  instrumentation. 

(a)  Autodata  9  machine  and  amplifier. 

(b)  Multichannel  pyrometer. 
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(7)  Turn  on  the  switch  to  the  feed  pump  to  recirculate 
water  in  feedwater  tank. 

(8)  Insure  nitrogen  level  in  cold  trap  is  at  full 
mark  (if  installed)  and  start  vacuum  pump. 

(9)  After  feedwater  tank  has  reached  a  temperature 
of  140  F,  insure  water  level  in  boiler  is  above 
low  level  mark  and  energize  boiler. 

(10)   Open  valve  DS-1 . 

b.   House  Steam  Operation 

Follow  steps  (1)  through  (4),  (6)  and  (8)  as  outlined 
above  for  the  boiler. 

2 .   Operation 

a.   Cooling  Water  System 

(1)  Open  valve  CW-1;  then  open  valve  CW-2  one  turn 
to  prime  the  cooling  water  pump,  keeping  valves 
CW-3  and  CW-4  closed. 

(2)  Energize  pump,  and  close  valve  CW-2.   Open  valve 
CW-3  one  turn  until  flow  is  established,  then 
open  valve  CW-4  to  purge  air. 

(3)  Open  valves  CW-3  and  CW-4  to  obtain  desired  flow 
rates . 

(4)  Vent  both  sides  of  the  12-ft  manometer. 

(5)  When  using  the  house  water  supply  remove  plug 
from  sump  and  open  valve  CW-2  with  valve  CW-1 
closed.   Follow  step  3. 
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b.  Steam  System 

(1 )  Boi 1 er  Operation 

(a)  When  boiler  has  reached  the  desired  pressure 
(approximately  3  psig)  open  valve  MS-1. 

(b)  Insure  valves  MS-6  and  MS-5  are  open. 

(c)  Open  valve  MS  -  3  to  obtain  desired  steam  flow 
rate  to  test  condenser.   Open  valve  MS-4  as 
necessary  to  maintain  boiler  pressure  at 
desired  level  . 

(2)  House  Steam 

(a)  Insure  valve  MS-1  is  closed.   Open  valve  MS-2. 

(b)  Follow  steps  (b)  and  (c)  for  boiler  use. 

c.  Condensate  and  Feedwater  System 
( 1 )  Using  Boiler 

(a)  To  collect  drains  in  test  condenser  hotwell 
operate  with  valve  C-l  closed.   After  test 
run  has  been  completed,  open  valve  and  con- 
densate will  drain  into  secondary  condenser. 

(b)  The  condensate  pump  is  operated  intermittently, 
when  level  in  secondary  condenser  dictates. 
When  pump  is  secured,  keep  valve  C-2  closed. 
When  pump  is  required,  start  pump  and  then 
open  valve  C-2.   In  this  mode  keep  valve  C-3 
closed. 
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(c)  While  feed  pump  is  running  (continuous  oper- 
ation) valve  FW-1  must  be  fully  open  and 
valve  FW-2  must  be  throttled  so  that  a  posi- 
tive flow  is  insured.   Valve  FW-3  is  a  sole- 
noid valve  which  is  actuated  by  the  boiler 
control s . 

(d)  When  boiler  is  energized,  valve  FW-4  must  be 
fully  open. 

(e)  Make-up  is  added  to  the  system  through  the 
top  of  the  feedwater  tank  by  removing  anode. 

(2)  Using  House  Steam 

(a)  Follow  step  (a)  for  using  boiler. 

(b)  To  pump  condensate  from  secondary  condenser 
hotwell,  start  pump,  and  open  valve  C-3.   In 
this  mode  keep  valve  C-2  closed. 

(c)  Delete  steps  (c)  through  (e)  for  using  boiler. 

Securing  System 
a.   Using  Boiler 

(1)  Close  valves  MS-3  and  MS-4.   Secure  power  to  boiler 
and  then  cl ose  MS- 1 . 

(2)  Pump  condensate  from  secondary  condenser  hotwell  to 
feedwater  tank.   Secure  valve  C-2. 

(3)  Secure  power  to  heater. 

(4)  Secure  vacuum  pump. 

(5)  Secure  cooling  water  pump  or  close  valve  CW-2  when 
using  house  water  supply.   Close  valves  CW-3  and 
CW-4. 
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(6)  Secure  instrumentation. 

(7)  Bottom  blow  boiler  to  remove  deposits.   Repeat 
twice,  blowing  from  high  water  mark  to  low  water 
mark. 

(8)  Secure  power  to  feed  pump. 

(9)  De-energize  individual  circuit  breakers. 
(10)   De-energize  circuit  breaker  on  control  panel; 

depress  stop.   Turn  key  switch  off. 

b.   Using  House  Steam 

(1  )   Close  valve  MS-2. 

(2)   Pump  condensate  into  return  line;  close  valve  C-3 
C3)   Follow  steps(4)  through  (6),  (9)  and  (10)  as 
outlined  for  procedure  using  boiler. 

4.   Secondary  Systems 

a.  Vacuum  System 

Vacuum  is  established  by  mechanical  vacuum  pump  and  is 
controlled  by  a  vacuum  regulator  mounted  on  instrument 
board  mounted  by  test  condenser. 

b.  Desuperheater 

Valve  DS-1  controls  flow  of  feedwater  (140°F)  to  spray 
nozzles.   Optimum  flow  level  is  between  15  and  20  per- 
cent flow  on  rotameter.   Condensate  is  collected  in  a 
small  tank  below  desuperheater  so  the  mass  flow  rate 
can  be  determined. 
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5.   Safety  Devices 

a.  Emergency  Power  Shut-Off 

To  secure  all  power  to  the  system  in  an  emergency, 
depress  the  red  button  on  the  right  side  of  the  main 
control  panel  . 

b.  Boiler 

(1)  The  mercury  switches  mounted  on  the  main  control 
panel  secure  power  to  the  heating  elements  of  the 
boiler  when  the  steam  pressure  exceeds  25  psig. 
Power  is  restored  to  the  heating  elements  when  the 
pressure  drops  to  approximately  15  psig. 

(2)  A  low  water  level  limit  switch  is  contained  within 
the  boiler,  and  when  the  water  level  inside  the 
boiler  drops  below  a  preset  level,  power  is  se- 
cured to  the  boiler  and  will  not  be  restored 
until  the  water  level  is  above  this  preset  height. 

(3)  The  relief  valve  mounted  on  the  boiler  is  set  to 
lift  at  30  psig. 
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APPENDIX  D 
SAMPLE  CALCULATIONS 

A  sample  calculation  is  performed  here  to  illustrate  how 
the  data  reduction  program,  Appendix  H  progresses  to  the  re- 
sults.  Tube  number  1,  run  number  2  at  60  percent  flow  on  the 
large  rotameter  was  selected  at  random  to  perform  this  analysis. 
This  tube  and  run  number  are  the  same  as  that  used  for  the 
error  analysis  in  Appendix  E. 

Section  2  of  this  appendix  corresponds  to  the  calculations 
performed  for  N  =  1  in  the  data  reduction  program.   The  calcu- 
lations for  N  =  2  will  be  demonstrated  in  section  3  of  this 
appendix.   The  water  property  calculations  are  shown  in  section  1 


INPUT  PARAMETERS 

Tube  Number 

Run  Number 

Tube  Inside  Diameter, 
Plain  End  (D.) 

Tube  Inside  Diameter, 
Test  Section  (D.) 

Tube  Outside  Diameter,  (D) 

Inside  Wetted  Perimeter  (Pw.) 

Outside  Wetted  Perimeter  (PwQ) 


1 
2 

0.0436667  ft 

0.02408333  ft 

0.0525  ft 
0.155833  ft 
0.183333  ft 


Test  Section  Cross  Sectional  Area(Ac)  0.0009375  ft' 

2 


Outside  Nominal  Surface  Area  (A  ) 
Wall  Resistance  (R  ) 


0.490874  ft' 


4.3761  x  10 


-5   hr-ft 


2  o, 


BTU 
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Cool ing  Water  In  (Tc  .  ) 

Cooling  Water  Out  (Tc  ) 

o 

Average  Cooling  Water 

Temperature  (T.  ,  T,  x 
v  v  b'   br) 

Steam  Vapor  Temperature  (T  ) 

Tube  Wall  Temperature  (T  ) 

Tube  Pressure  Drop  (AP  ) 

m 

%    Flow 

Tube  Inlet  Contraction  Factor 

Tube  Outlet  Expansion  Factor 

*See  Appendix  I . 


71 .24°F 


80.74UF 


75.99°F,  533.59°R 


151.07UF 
554°R 

19.56929  psi 
60 

Kc  ) 

>    0.33 
Ke 


Section  1,  Water  Properties 


yH20 


H„0 


exp[( 0.004606532) (533. 59)+ (4759. 5941 )/( 533 . 59) 

-  10.59252566] 


2.19333 


lbm 
ft-hr 


0.32159931  +  (0.000697989)(75.99)-(0. 12506  x  1 0"5  )  ( 75 . 99 )2 


-(0.2072  x  10"10)(75.99)3 


0.36741 


BTU  _ 
hr-ft-°F 


62.707172  -  (0.0043955304)(75.99) 


62.10706 


-  (0.000046076921)(75.99) 

lbm 
ft3 
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1.0121559  -  (0.00024618473)(75.99) 

+  (0.10282155  x  10'5)(75.99)2 


0.99939 


BTU 
lbm°F 


GPM  x  p   =   (ll.28)(60)(^g)  (62.10700) 


5647.22266 


Ibm 
hr 


1  .57 


Ibm 
sec 


Pr 


uc, 


(2.1 933  )  (0. 99939  ) 
.36741 


5.966 


Section  2,  Plain-End-Tube  Reduction 


1.   Determination  of  cooling  water  velocity 


4m 


PttD. 


TS 


pAc 


(4)(5647. 22266) 


5647.22266 


(62. 10706)  (tt)  (  .  043667) 


(0.0009375) (62. 10706) 


60716.023 


ft 
hr 


TS 


96989.04 


11 
hr 


16.87 


ft 
sec 


26.94 


ft 
sec 
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2.   Determination  of  Mass  Flow  Rate  per  Unit  Area 


4m 


ttD 


pv 


(62. 10706) (60716 .023) 


3,770,893.659  -l|m 


ft  -hr 


1047.47 


Ibm 


ft  -sec 


3.   Determination  of  Reynolds  Number 


Re  = 


Ho0 


(0.0436667)(3770893.659) 
2.1933 


Re  =   75,075.22 


4.   Determination  of  Overall  Heat  Transfer  Coefficient 


m  c       T  -  Tc. 


(5647. 22266) (0. 99939) 

0.490874 


In   [ 


151  .07  -  71  .24-, 
151 .07  -  80.74J 


1456.83 


BTU 
hr-ft2  °f 
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5.   Determination  of  Corrected  Overall  Heat  Transfer  Coefficient 


w 


1 


1 


1456.83 
1556.03 


-  4.3761  x  10 

BTU 
hr-ft2  °F 


6.   Determination  of  Friction  Factor 
0.046 


Re 


0.2 


0.046 


75075.22 


f   =   0.00487 


0.2 


4feG2(nL) 


AP 


AP 


AP 


P  2gc 

(4)  (0.00487)  (3770893.659  )2  (-^ll^-) 
(62. 10706) (2) (4. 1 7  x  1 08)  ( 1 44) 


1  .072  psi 


AP 


P  V 


TS 


exp/con 


2g 


[Kc  +  Ke] 


(62.107172H96989.04p      [0.33] 
(2)  (4. 1 7    x    108)    (144) 


AP 


exp/con 


1 .6054    psi 
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AP 


TS 


TS 


TS 


AP  -  AP   -  AP    , 

s     exp/con 


19.56929  -  1 .072  -  1 .6054 
16.89  psi 

pAPTs2gc 

4G2  f|Jr) 

(62.107172)(16.89)(2)(4.17  x  108)(144) 
(4)(3770893.657)2  (^ggj) 

0.032246 


7.   Determination  of  Wilson  Plot  Parameters 


(a)   Ordinate 


U 


1 


1456.83 
(b)   Abscissa 


1 Uu  n 

0.8  "  1/3  /   M2Ux  0.14 


Reu.o  pr./,  ( 


y 


■) 


w 


w 


exp[ (0.004606532) (554)+ (4759. 5941/ (554) 

-  10.59252566] 
Ibm 


w 


1 .7347 


ft-hr 
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(75075. 22)°-8(5.966)1/3(ftl|f^)0-14 


6.71086  x  10 


-5 


8.   Determination  of  Sieder  Tate  Constant 

c   .  ^ 
Li     Mk 


1.76639,   from  linear  regression  subroutine 

0.0525 

(1  .  76639)  (0. 36741 ) 


c   =   0.080895 


9.   Determination  of  Inside  Heat  Transfer  Coefficient 


;ik  Re0.8  prl/3  (!V)0.14 


Di 


0.080895 


w 


0.8 


=  (q  ogggg^)  (0.36741)(75075.22)u'°(5.966) 

2.19333N0.14 


1/3 


( 


1 .7347 


-) 


hi 


10142.384 


BTU 


hr-ft2  °F 


10.   Determination  of  Outside  Heat  Transfer  Coefficient 

1 

Un  "  Rw  "D.h. 


1 


6.86422  x  10'4-  8.75223  x  1  0"5-  (Q  ■  04^667  )  ( 1  01  42  .  384) 


2081 .777 


BTU 


hr-ft2  °F 
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11.   Determination  of  Nusselt  Number 


Nu 


h.D. 


1205.43 


12.   Determination  of  Stanton  Number 


St 


Nu 


RePr 


1205.43 
(75075. 22) (5. 966) 


2.691  x  10' 


13.   Determination  of  Performance  Factor 


TPF 


2xJ 
fTS 


St  Pr 


2/3 


8.8529  x  10' 


TPF 


(2)(8.8529  x  10"3) 
.032246 


TPF   =   0.5491 


Section  3.   Tube  Reduction  Based  Upon  Hydraulic  Diameter 


1.   Determination  of  Velocity 


TS 


26.94  ft/sec,  96989.04  ft/hr 


*  See  Section  2  of  this  Appendix 
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2.   Determination  of  Mass  Rate  Flow  per  Unit  Area 


G   =   Ac" 


5647.22266 
.0009375 


6023704.171 


lbm 


ft  -hr 


1673.25 


lbm 


ft  -hr 


3.   Determination  of  Reynolds  Number 


Re  = 


DhG 
yH„0 


(0.02408333)(6023704.171 ) 
2.1933 


66,142.74 


4.   Determination  of  Overall  Heat  Transfer  Coefficient 


U   =   1456.83  ^ — 

n  hr-ft2  °F 


5.   Determination  of  Corrected  Overall  Heat  Transfer  Coefficient 


U   =   1669.73  BT~  - 

C  ,        ,,Z    Or 

hr-ft   F 


6.   Determination  of  Friction  Factor 


AP 


TS   =   16.89  psi 


*See  Section  2  of  this  Appendix. 
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TS 


pAPTs2gc 

uh 


(62. 107172)(16. 89)(2)(4. 17  x  108)(144) 
(4)(6023704.171)2  (QtQH^m) 


TS 


0.006968 


7.   Determination  of  Wilson  Plot  Parameters 

(a)   Ordinate 

Y  =   4- 

Un 

Y  =   6.86422  x  10"4 


(b)   Abscissa 
X   = 


1 


Re0.8  Rrl/3(V  0.14 


8. 


,,,,,  7/1O.8  "  c  nccl/3  ~  ,2.  1933N0.14 
66142.74    x  5.966     x  U — jttj) 

X   =   7.4266  x  10"5 

Determination  of  Sieder  Tate  Constant 
AnD, 


ci 


Pw.jL-j.5Mkb 


(0.490874)(0. 0240833) 

(0.155833)(3)(1.59616)(0. 36741) 


c.   =   0.04312 
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9.   Determination  of  Inside  Heat  Transfer  Coefficient 


c,k 


0.8  D  1/3  ,  H2°s0.14 


1_  Reu-a  Pr,/3  (-*-) 


h 


P. 


(0.04312 ) (0. 36741 ) 
(0.0240833) 


(66142. 74)°-8(5.966)1/3 


,2.1933x0.14 
M  .  7347 ; 


8857.72  hrBTfUt,0p 


10.   Determination  of  Outside  Heat  Transfer  Coefficient 


Pwrt    L  Pvr 

o  o 


Pw 


AnUn  Pbar      w        Pwihi 


1 


(0-18333H3)    6>86422    x    1Q-4    _    0.1833         8>75223x10-5 
(0.490874  0.1695815 


1  .12896    x    10"4 


155833 


h        =      1846.2   5!" 

0  hr-ft2    °F 


11.      Determination    of   Nusselt    Numbe. 


Nu 


hiDh 


(8857. 72)(0. 0240833) 
0.36741 


Nu 


580.61 
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12.   Determination  of  Stanton  Number 


St  = 


Nu 


RePr 


St 


580.61 

(66142.  74)(5. 966) 


St 


1 .4714  x  10 


13.   Determination  of  Tube  Performance  Factor 


TPF  = 


2  x  J 
fTS 


St  Pr 


2/3 


1 .4714  x  10"3  5.9662/3 


4.84  x  10 


-3 


TPF  = 


(2)(4.84  x  10~J) 
0.006968 


TPF  =   1 .39 
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APPENDIX  E 
ERROR  ANALYSIS 

The  basic  equations  used  in  this  section  are    reproduced 
from  Pence  [13],   The  general  form  of  the  Kline  and  McClintock 
[26]  "second  order"   equation  is  used  to  compute  the  probable 
error  in  the  results.   For  some  resultant,  R,  which  is  a  func- 
tion of  primary  variables  X] ,  X«,  ...  X  ,  the  probable  error 
in   R,   <5R  is  given  by: 


6R 


^it  5Xi>2  +  <4k  6X2>2  + 


1 


+  (jf   ^n)2]V2  (E-l) 


where  6X -j  ,  6X2,  ...,  6Xn  is  the  possible  error  in  each  of  the 
measured  variables. 

The  overall  heat  transfer  coefficient  is  given  by 
equation  (4),  in  Chapter  III  as: 


T   -  Tc . 
,i    _   m  cp   ,    r  v i-, 

Un   '  -nT      ln   [TV  -  Tc^   • 


(4) 


By  applying  equation  (E-l)  to  equation  (4),  the  following 
equation  results: 


5U. 


[(^r)2  +  (4^)2  +  (^)2  +  ( 


6Tv  (TcrTco) 


An 


cp 


m 


(Tv-Tci>(VTco)ln— 


VTco 


6Tc  2 

< °-T^TT)       } 

(Tv-Tco)lnT^TEl 


1/2 


(E-2) 


182 


The  following  are  the  values  assigned  to  the  variables 


Sep  =  0.001  BTU/lbm  °F 

6mc  =  0.01  m  Ibm/hr 

6TS  =  0.2  °F 

6Tc.  =  0.2  °F 


6Tc 


0.2  °F 


6An 


0.008  ft' 


The  probable  error  in  the  Reynolds'  number  is  given  by: 


5  Re 
Re 


[(6G)2  +   ^2  +  C^l)2]1/2    t 


(E-3) 


where , 


6G 
G 


[(04Ji,Z  +   (^i,2]l/2 

m  i 


(E-4) 


6y 
<SD 


l 


0.01  lbm/ft-hr,   and 


0.00042  ft 


The  probable  error  in  the  coefficient  c.  is  given  by 


4c, 


SD 


[(^)2  +  (4^)2  +  (^)2:1/2     . 


slope 


(E-5) 


where 


<5D. 


0.00092  ft 


6k  =  0.001  hrB^.oF   ,  and 


6  slope   =   0.065  slope 
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The  probable  error  in  the  inside  heat  transfer  coefficient 
is  given  by: 


6hi   .   [C^}2  +  (^1)2  +  (^i)2  +  (0_^5Re}2  +  (0.3336Pr}2 


hi 


Re 


Pr 


u  /U. 


w 


where ; 


5k 
6D 


0.001 


BTU 


6CT 


6Re 


hr-ft-  F 
0.00042  ft   , 

from  equation  (E-5)  , 

from  equation  (E-3)  , 


6Pr 


0.01   ,   and 


V 


0.01 


The  probable  error  in  the  outside  heat  transfer  coefficient 
is  given  by: 


6h 


C  ( 


611 


0 


)2+  ( 


5R 


Dihi 


(n1-  -  R 


w 


Dihi 


+  ( 


D     6h. 


1/2 


J--  R 
%.    w 


D 


Dihi 


(E-7) 
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where: 

6U 

<5R 


w 


6h 


from  equation  (E-2)   , 

8.75  x  10~6  hr-ft2-°F/BTU,   and 

from  equation  (E-6)  . 


The  probable  error  in  the  Nusselt  number  is  given  by 


6Nu 
Nu 


[(;v  ♦  i64f)2  ♦  <^)2  ]i/2  -, 


(E-8) 


where 


6h. 


as  found  in  equation  (E-6)   , 


<5D 


6k 


0.00042   ft,    and 
BTU 


0.001 


hr-ft-F 


The  probable  error  in  the  Stanton  number  is  given  by: 


6Sfr 
St 


r  (&M\2        (6*e.\z        (&lL\2   -i1/2 
L  K   Nu  ;     l"Re  ;     K   Pr  ;   J 


(E-9) 


where: 

5Nu 
5Re 
6Pr 


as  found  in  equation  (E-8), 

as  found  in  equation  (E-3),  and 

0.01   . 


The  probable  error  in  the  friction  factor  is  given  by 


6f 

f 


6AP 


TSx2 


6Gx2 


6Lx2     ,6pv2 


C  (TP=r)£  +  (2  TT>  +  (l-}  ♦  (  p  ' 


♦  (-DV  ]1/2 

uh 


(E-10) 
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where: 

SAP 

5G 
6L 

<5f 


TS 


0.009  AP 


as  found  in  equation  (E-4), 


0042  ft   , 


0.05  Ibm/ft  ,  and 


6DL    =   0.00042  ft 


The  probable  error  in  the  Tube  Performance  Factor  is 
given  by: 


6TPF 
TPF 


r  c^)2  +  4^)2  +  (^)2  ]1/2,  (E-iD 


where 


St 
Pr 
f 


as  ijund  in  equation  (E-9), 

0.01  ,  and 

as  found  in  equation  (E-10) 
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For  Run  #2  at  60  Percent  Flow 


6^   _  ,.008       ,2       .  .2  ,.01nu2        (  ,    U)(-9  45)  . 

Un    '  \490874;         uuuu  l_F~J  (79.83)(70.335)ln(1.135)j 

n9  2     1/2 

v  (70.335)ln(1.135)y      J 


^-n  =         0.030 
Un 


Un,   60%       =       1456    '     50       hf^op 


6G     =         r     012  +   (-00042    .2  -,1/2         =     n  Q2 
G  L   '  \024083;     J  ^-^ 


6Re 
Re 

6Re 
Re 


r  (  Q2)2  +  (.^P-L-)2  +  c00042  )2  i1/2 

L   l-^l     +   v2  -,g;     +   ^.024083;     J 


027  .  *  .    Re      =      66142    ±    1786 


5ci  r    f.00042v2         ,   n«i2  ,.001      %   -,1/2 

"c-  [   (^525_)     +     ('065)  (^674T)   J 


6ci 

— -  =       .066 

ci 


ci,60%     "     '°431   "     -°028 


6h 


,    .01    .2      .      f0.14  x   .01,2   -,1/2 
+  (5T966j  [       1.034     J     J 
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-r-1       =  .068 


hi ,  60  %     =     8861   "  603  EFT^F 


6ho  C"(  50 


h°  167Q2(W0-8-75x1Q'5-886i°x250437)2 


+  /  8.75  x  TO"6  .  2 

^T670  "  8'75  x  10"J  "  8861  x  .0437^ 


^0525__     x  _068  2    1/2 


8861  x  .0437 


,   /   OOP  I  X  .UHJ/ x      -i 

1   -  8.73  x  ID'5     °'525 


1670   «.,.,  a  ,M     8861  x  >Q437 


-r-°-    =     0.058 
ho 


BTU 


ho,  60%     "     1846  "   107  hr-ft2  of 


<5_Nu 
Nu 

Nu 


r    t   n^\2       ^-00042   ,2       ,    .001    ^2   -,   1/2 
Nu  [   ('068)     +  (^24083)     +  (3674T)     ] 


=     0.070 


Nu       =     581   -  41 
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&■    =     [   (.07)2+  (.027)2+  (^Ig)2  ]1/2 


^    =       075 
St  'u/0 


St       =     .001472  -   .000110 


f      =     [   (.009)2  +   (2  x   .02)2  +  (^-)2  +   (^§p)2 


,.00042   ^2  -,   1/2 
\  024083 '     J 


^      =     .045 


00697  -   .00031 


<5TPF  r    ,   n7,>2        ,4       .01 ^2       ,   n/lc,2  n  ]/2 

TPT    "     [   (-075)     +  (3  x  5T966)     +  ('045)     ] 


6TPF 
TPF 

=     0.087 

TPF 

=     1.39  - 

.12 
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APPENDIX  F 
AREA  RATIOS,  SAMPLE  CALCULATIONS  AND  RESULTS 

To  determine  the  relative  performance  of  the  General 
Atomic  tubes  as  compared  to  smooth  tubes,  it  is  necessary  to 
define  a  method  for  comparing  the  augmented  tubes  to  the 
smooth  tube.   Such  a  comparison  is  found  in  Bergles  [3]  and  is 
described  in  Chapter  III  of  this  report  for  the  case  where 
Roy+  =  0.   By  solving  the  equations  as  shown  in  Chapter  III, 
it  is  found  that 


Re 


=  {«<~f 


2  Nu  /Pr 
a 


0.4 


(F-l) 


In  this  case,  for  a  given  cooling  water  flow  rate,  Re  .  f  , 

a   a 

Nu3  and  Pr  are  all  found  in  Tables  14  through  19  depending  on 
a 

which  tube  is  being  investigated.   Also,  knowing  Re   permits 
the  solving  of  the  area  ratio: 


046  Re 


2.8 


Rea  fa 


(F-2) 


For  run  2,  45  ,  the  summary  of  area  ratios  are  listed  below 


Rea 

fa 

Nu 
a 

Pr 

Re 
s 

VAs 

40176 

.03  5 

719 

5.66 

56191 

.40 

64696 

.033 

1058 

5.79 

92332 

.41 

89868 

.032 

1376 

5.84 

130750 

.42 

94931 

.031 

1437 

5.84 

136720 

.42 
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For  run  2,  45  ,  the  summary  of  area  ratios  are  listed  below 


Rea 

fa 

Nua 
a 

Pr 

Res 

A  /A 
a   s 

25872 

.038 

522 

5.66 

35510 

.39 

37416 

.032 

715 

5.90 

48827 

.37 

61178 

.029 

1067 

6.03 

79902 

.37 

83894 

.027 

1387 

6.17 

109089 

.37 

95463 

.026 

1541 

6.20 

123396 

.37 

For  run  9,  60  ,  the  summary  of  area  ratios  are  listed  below 


Rea 

fa 

Nu3 
a 

Pr 

Res 

VAs 

26384 

.064 

531 

5.53 

46836 

.47 

38331 

.059 

725 

5.74 

67897 

.47 

61813 

.059 

1080 

5.96 

114781 

.48 

79970 

.062 

1332 

5.99 

156065 

.50 

By  expanding  on  Bergles1  [3]  procedure  to  include  the  case 
where  Rovf    t    0,  as  shown  in  Chapter  III  of  this  report,  the 
following  equation  can  be  found. 

fa  yl    c    pDi  i/q    ] 
v.   =   [  a   a     (-r1)'75]  2"3  (F-3) 

s      U   .046    p 
a 

again,  for  a  selected  v  .  values  of  v  .  f  ,  and  U   can  be  found 

a  s    a       a 

Knowing  v  , 


U 


c  V^ 


(F-4) 


where  c   =   251  is  a  constant  and  was  found  as  shown  in 

Chapter  III. 

An  area  ratio  is  now  found  from 
A.     II 


(F-5) 
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For  run  8,  30  ,  and  using  the  procedure  outlined  in  Chapter 
III: 

D.    =   0.53  in 

Vavg    =   2.17  lbm/ffhr 

p    =   62.1  lbm/ft3 
Tci    =   65°F 

By  applying  the  temperature  correction  factor  to  U   as  found  in 

reference  [21],  U,  is  found  from: 

a 


.9T 


(F-6) 


For  v 


2.73  fps,  substituting  the  required  values  above  into 


equation  (F-3),  it  is  shown  that 


r,(.053)(2.73)3  251  v   ,(62.1  )  (  .  53  )  (3600  ) ]  7  5  ^ttU- 
l[         (1050)  (  .046) — }   ( (2.17)(12) ]  J  2  "  3 


4.39  fps 


By  substituting  this  value  into  equation  (F-4), 


U    =   251  V4739 
s 


526 


BTU 


hr-ft2  °F 


Using  the  results  of  equations  (F-4)  and  (F-6),  the  area  ratio 
is  found, 


526 
1050 


0.50 
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A  Reynolds  number  corresponding  to  v   is  found  from 


Re. 


p  D  •  v 
H  l   s 


.53 


(62.1)(-1-2~)(4.39)(3600) 


2.17 


Re 


19975 


are 


For  run  8,  30  ,  the  remaining  area  ratios  for  Rov+.  f    0 


a 

f  ps 


s 

f  ps 


Re 


5.47 

8.63 

.58 

39268 

8.20 

13.22 

.67 

60153 

13.67 

23.31 

.79 

106065 

19.14 

34.36 

.91 

156344 

21  .87 

39.35 

.93 

179050 

Similarly  for  run  2,  45°,  the  summary  of  area  ratios  for 


ext 


f    0  are  found  to  be: 

Tc,  =  70' 


a    n 


va 
f  ps 

vs 
fps 

Aa 
As 

Res 

2.90 

4.52 

.54 

21012 

8.59 

14.95 

.74 

69497 

14.12 

26.81 

.92 

124630 

19.75 

39.88 

1  .05 

185387 

20.86 

42.32 

1  .07 

196730 
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and  for  run  9,  60  ,  the  summary   of  the  area  ratios  for 
R   .  f   0  are  found  to  be: 


a 
fps 


s 
fps 


Re. 


2.73 

4.69 

.51 

22157 

5.47 

11.14 

.69 

52629 

8.20 

17.24 

.76 

81448 

13.67 

32.08 

.92 

151557 

17.77 

44.97 

1.03 

212454 
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APPENDIX  G 
CALCULATION  OF  WALL  THICKNESS  AND  WALL  RESISTANCE 


There  were  three  methods  considered  for  determining  the  wall 
thickness  which  are: 

1 .   Anul ar  Area  Method 

A  nominal  outside  area  is  found  using  an  outside 
diameter  of  5/8  inch  and  substituting  into  equation  (G-l  ) 
written  as: 


D 


i  -fo 


2   4  Amet 


(G-l) 


where  Amet  is  found  in  Table  3.   The  wall  thickness  is  found  by 
substituting  D.  and  D   into  equation  (G-2),  written  as: 


D-D. 
o    i 


(G-2) 


2.   Displaced  Volume  Method 

The  volume  despl acemen t ,  V~,  of  the  tube  is  measured. 

An  average  perimeter,  P,  is  determined  from  equation  (G-3)  below 

After  substitution  of  V^  and  F,  the  wall  thickness  is  found  from 

equation  (G-4)  below: 

P  +  P. 
P    =    °  o       3      ,  (G-3) 


V, 


P  L 


(G-4) 


3.   Indirect  Measurement  Method 

A  photograph  of  the  cross  section  of  the  experimental 
tube  is  enlarged.   The  wall  thickness  could  generally  be 
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measured  from  the  enlarged  photograph.   However,  as  seen  in 
Figure  13,  the  wall  thickness  varies  around  the  periphery  of 
the  tube. 

Method  1  was  selected.   Method  2  was  rejected  because  it 
was  felt  that  the  displaced  volume  would  be  extremely  difficult 
to  measure.   Method  3  was  rejected  because  of  the  non-continuous 
variations  of  the  wall  thickness. 

The  wall  thickness  for  each  tube  is: 
tube  1 

=   0.4811  in 


Di 


w 


tube  2 


tube  3 


w 


Di 


w 


0.07195  in, 

0.5006  in 
0.062215  in,  and 

0.4942  in 
0.0654  in. 


The  wall  resistance  is  calculated  by  solving  equation  (G-5), 
written  as: 


R 


w 


t 
w 

k    ' 
W 


(G-5) 


where  K   is  found  to  be  137  BTU/  hr-  F-ft  in  the  Handbook  of 
Tables  for  Applied  Engineering  Science  [27]. 
For  tube  number  1 : 


R   =   4.3761  x  10 
w 


-5  hr-ft2-°F 
BTU 
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For  tube  number  2 


w 


3.7844  x  10 


5   hr  ft2  °F 
BTU 


For  tube  number  3 


R   =   3.9781  x  10 
w 


-5   hr  ft2  °F 
BTU 
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APPENDIX  H 

DATA  REDUCTION  PP.OGP.AM 


IV  G  LEVEL   21 


MAIN 


CATS  =  78075 


C3M3/C3 


C  DEFINITION  OF  SYMBOLS  US£D  IN  THS  PROGRAM 

f 

C  AC=TES7  7L3E  CROSS  SECTIONAL  AREA  IN  SQUARE  F6E7 

C 

C    ACTL=A    CHANGING    VARIABLE 

r 

C    ACT2=A    CHANGING    VARIABLE 
r 

C    ACT2=A    CHANGING    VARIABLE 

C 

C    ACT4=A    CHANGING    YARIABLF 

C 

C  ANaARSA  Or  STEAM  SICE  OF  TES7  TUBS  IN  SQUARE  FEET  3ASED  CN  5/8  IN  CD 

C  ?.--Y  INTERCEPT  AND  SLOPE  0=  EQUATION  OBTAINED  SY  LINEAR  REGRESSION 
r 

C  CiaSIECSP  TATE  CONSTANT 

C  CP-CCCLING  WATER  SPECIFIC  FEAT 

C  CCRR=THE  NGN  RECOVERABLE  LOSS  COEFFICIENT;  DUE  TO  EXPANSION  AND  CONTRACTION  IN 
C  THE  TEST  SECTION  OF  THE  TEST  TUBE 

r 

C    DELY=ThE    CUTPUT    ARRAY    OF    THE    DIFFERENCES    BETWEEN    Y    AND    YO 

C 

C    DELP  =  'EST    TUBE    PRESSURE    D?.CP 

f 

C    0ELPC=PlAIi\    ENO    PRESSURE    CROP    CORRECTION 

r 

C    CELP=C=NON    RECOVERABLE    PRESSURE    DROP    CORRECTION 

C 

C    DI=INSIOE    DIAMETER    CF    TUEE    IN    FEET 

r 

C    DH=HYCRAULIC    DIAMETER   OF    TUBE    IN    FEET 

C    OLTP=?RESSURE    DROP    ACROSS    TEST    SECTION    OF    TEST    VU3E 

C    CU  =  OuTSIDE    DIAMETER    OF    TUEE    IN    FEET 

C 

C    CCTHsCGGLING    WATER    PASS     FLOW    RATE 

C    =«PERC3NT    CF    FULL    SCALE    ON   TH  =    ROTAMETER    OF    WATER    FLCblNC    THROUGH    THE    TU3E 

C    FC=CALI3RATICN    CORRECTION    FOR    COOLING    WATER    FLOW    THFOUCH    ROTAfETER 

C 

C  FRIFAC=TEST  SECTION  FRICTION  FACTOR 

C  FSMTU-^?LAIN  END  FRICTION  FACTORC 

C 

C    C-  =  ,'ASS    FLCW    RATE    PER   UNIT    AREA 

C 

Z    HI=INSIOE  HEAT  TRANSFER  COEFFICIENT 

C  HO=OUTSICE  HEAT  TRANSFER  COEFFICIENT 

C 

C    H2G=C00LING   WATER    CCNOUCTIVITY 

0    L5C?L2  =  LINEAP.    REGRESSION    SUBROUTINE  (SEE    REFERENCE  ) 

C     PS  A  R=  AVE  RAGE    OF    PO    ANO    P'.J 

C    PNU=NUSSELT    NUMBER 

C 

C     PO=OUTSIDE    WETTED    PERIMETER 

C 

C    PR=CGCLING    WATER    PRAMOTL    NUMBER 

C.    PW=WETTED    PERIMETER    OF    COOLING  rfATER    SIDE    OF    TEST    TL3E    IN    FEET 

C    RATMHL*RATI'J    OF    WATER    VISCCSITY    TO    TU3E    WALL    VISCOSITY 

C 

C    PE=COCLING    WATER    REYNOLDS    NUMEER 
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C    REAN=TU3E    PERFORMANCE    FACTOR 

C    RE2=PLAIN    END    REYNOLDS    NUMBER 

C 

C    RE3=TEST    SECTION    SEYNCLOS    NUMEER 

C    RHC=COCLlNG    WATER    DENSITY 

C 

C  RHOP=CCJLING  WATER  DENSITY 

C 

C    RW=W£LL    RESISTANCE 

r 

C    S3=ESTIMA7E0    ERROR    CF    B 

C 

C     ST=STANTCN    NUMBER 

t    Tc=AVERAG£    OF    INLET    AND    OUTLET    COOLING    HATER    TEMPERATURES 

C 

C    T3R=T9    IN    CEGREES    RCENGTEN 

C    TCI=CCGLING    WATER    INLET    TEVPE  RATURE    IN    CEGREES    F 

r 

C  TCC=COCLING    WATER    CUTLET    TEMPERATURE    IN    CEGREES    F 

C 

C  TS=VAPGR    TEMPERATURE    MEASURED    AT    INLET    TO    THE   CONDENSER 

C 

C  T!*  =  TEST  TL3E  WALL  TEMPERATURE  IN  DEGREES  RCENGTEN 

C 

C  UC=CORRECTEC    OVERALL    KEAT    TRANSFER    COEFFICIENT 

C    UN=QVERALL    HEAT    TRANSFER    COEFFICIENT 

C 

C    V=COCLING    WATER    VELOCITY 

C 

C    VTS  =  VELOCITY    CF    CCCLING    WATER    I.J    TEST    SECTION    OF    TEST    TUBE 

C    WALMHU=VlSCOSITY    AT    TEST    TUBE    WALL 

r 

C    W*HU=CCOL;NG    WATER    VISCOSITY 

C 

C  Wl*VARIABLE  WEIGHTS  USEO  IN  LINEAR  REGRESSION  SU3R01TI.ME 

C 

C    X=MA7RIX    CF    COORDINATE    POINTS    PASSED    TO    LINEAR    REGRESSION    SUBROUTINE 

C 

C  Y=MATRIX  CF  ORDINATE  POINTS  PASSED  TO  LINEAR  REGRESSION  .SUBROUTINE 

C 

C  VCGGRO=P.\'U/PR**  (1/3  )/'.  U/UW  )**  .14 

C 

C    YO=ThE    OUTPUT    ARRAY    OF    DEPENDENT    VARIABLES 

REAL    *    3    X.Y  .WI,  Y0.DELYT3,  S3 

D  I. ''ENSIGN  TS(20),7a(20).TBR(20)  .Tw(2J)f'JC(20),UN(20)  fRE(20)  ,F(20)» 
IFC120) ,V(20 ) ,PR(20) ,X(20) .Yt20) ,RATKHU(20) ,WI(20)»Y0(20),0ELY(20), 
23(2),S3(2)»TCI(20)»7C0(20i,Hl(20)  »hCM2C)  ,H20!<(20>  .0(20)  ,PNU(20)  .  FR 

3iFAC<20)  ,C\ (20)  ,YC00R0(2C  ),ST(20) ,FCHECK( 20) ,0ELP<20) .RE2(20) ,C3J( 
42 C), RE AN (20)  , RHG?(20)  ,VTS  (20)  ,  \zb  (  2  0  )  ,CLTP(  20  ) 
READ (5,1 J2) JK 
DO    ICC1    L  =  l  ,JK 
READ! 3  ,102) M 
102    FORMAT (113) 
C 
C    ZERO    ALL    MATRICES 


C 


DO    1    1=1,20 
TS( I )=0.0 
Ted  )*0.3 

T3.=  {  I  )=0.0 
TV.il »=0.0 
U\( I )=J.J 
UC(I )=0.0 
RE(I 1=0.0 
F(  I)=0.0 
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FCCI  )=0.0 
V(I)=C.O 
PP(I)=C.O 
X  (I)=0.0 
Y<  I)  =  C.O 
RATMHUl  I)«0.0 
w  I  ( !  1  ~  0  .  0 
Y  0  {  I  )  =  C  .  0  . 
0ELY( I  >=J.O 
TCid  1=0.0 
TCOU  )=0.0 
Hill ) =0.0 
HO (I ) =0.0 
H2C<(  I)=J.O 
G«  U  =0.0 
PMl  I  >  =  0.0 
FR  IF.'.c:  I)=J  .0 
C  K  i  I  )  -  0  .  0 
YCCGRCt  I)  =0.0 
ST(! )=0.0 
FOECM  15=0.0 

cel  = < : >=o.o 

RhCP(II=0.0 
V75(I 1=0.0 
RE3( I )=0.0 
DLT?< I) =0.0 
I    CCNTINUE 

am  =0.0 

3<2)=C.O 

SB(1)=0.0 

Sc<2)=0.0 
C 

C    REAO    IN    THE    RArf    DATA    AND    ALL    CONSTANTS 
C 

REflO(5i  l031AC,PV».00tPC»DH,ANtRW 

103  FGR'-'.AT(6F10.7tcl2.5) 

RE AD (5, 1041     ITS(K) ,T3(K)  ,T3R{K) ,TW(KJ»TC. .K) ,TCQ(K) ,F{K) ,FCIK)»K*1 
1 ,  V ) 
READ(5,104J     (DELP(K) ,K=1,MJ 

104  FORMAT  U-10.5) 

RE  A3  I 5,105101  ,CCRR 

105  FC=MAT(2F13.7J 
WR  I7=(6,29G) 

299    FGa-IATCl't  IX, * LNCORRECTEC    PRESSURES    AND    WATER    FRCP^RT  IES '// ) 
CC    10    1=1, M 
'/.SIT=(6,299)DELP(I  ) 
10   COM  i'^'JE 
290    FC."MAT(1X,F10.5) 
>i  =  C 

PEAR=(P0+PW)/2 
C 
C       CALCULATE    WATER    PROPERTIES 


C 


WPITE  (6,200)PW 

200  FCRKA7UX,el2.5  I 
9  0    N=N+i 

0C2    1*1, M 

IriN.RT.l  jr,Q    12    99 

ACU  =  J.  3  046  06  53  2"^  R(  I  ) 

AC?2=475S.5941/7QRl! ) 

W?-,hU=EXP(ACTl+ACT2-10.  59252566) 

WRITE(6,201 JWXHU 

201  FORMAT!  LX,F10.5J 

ACT1  =  C.0006V7989*T3(  I  ) 

AC72=c.i2  5:en-05jtTe  d  )*==2 

mC7?=C.£U7  2E-10*TQ( ! )**3 

H^Cn  (  I)=0.32159531+AC71-ACT2-ACT3 

WRITElo  f.iOl  )H2CK(I  ) 

ACTi=C.4j'.'i553'J4£-02*TB(  I  ) 

ACT2=:.46  07  69215-0'V--"r5(  !  )  **2 

RHG=62.707172-ACT1-ACT2 

WRITE(6,201)RHQ 

RI-CP(  I  l=RHO 
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ACTl=C.24fcl8^72E-03*TE(I ) 

AC72»C. 102821 55E-05*TB< 1 1**2 

CP  =  i. 3121559- ACT1+ACT2 

WRITE  15,201  )CP 

ACT1  =  C.C^46C6  53  2:ST^^.,  (  I  ) 

ACT2=4759.5941/TW< I ) 

WALMHU=5XP(ACT1+ACT2-10. 59252566) 

WRITE(6,201lhALPHU 

RAT,MHU(I)»IW:''HU/WALMHU)*-*3.14 

DC7M=F(I )*FC< I) *RH0*1 50. 7916667 

WRI7E(6t201 1D0TM 

V(  I)  =4. *D07M/<RK3*3.  141593*01**2) 

V7S(  I  )=DC)TM/(RHC*AC> 

RE  ( I  )=RK3*DI*V(  I J/WMHU 

AC71=i. 27324/01**2 

Gil) =ACT1*03TM 

RE2(I 1=01*G{I)/WMHU 

RE3I I)=D07M*CH/<AC*WMHU) 

PR  (I  )=CP-WMI-U/n20K(  I  ) 

c 

C    CALCULATE    THE    OVERALL    HEAT    TRANSFER    COEFFICIENTS 

C 

act:=tscj-tci(I) 

AC72  =  7Sm-7C0l  I) 

AC73=4CTl/iC72 

AC74=ALCGUCT3J 

UN  (I  l*DGT.M*CP*ACT4/AN 

AC71=1./UN( I) 

ACT2=«CT1-RW 

UC(I >=1./ACT2 

Y ( I)=ACT1 

GC    70    34 
C 

C  CALCULATE  THE  PRESSURE  DRCP  AMD  FRICTION  FACTOR 
C 

99  P.  =  <I) =  RS3< I) 

V  ( I )  =  V  T  S  <  I ) 

34  CONTINUE 

IFCRE21I1.LT. 3. E041G0  TO  36 
FSMTh=0.046/R52 (1 1**0.2 
GO  TO  3  5 
3  6  Fo^TH=3.079/RE2 (11**0.25 

35  ACTl  =  ^,.C-FSyTl-^G(I  1**2/01 
ACT2»RH0P  ( I  3  -:  1 .  2C096  E  +  ll 
D£LPC=1.260*ACT1/ACT2 
IFiN.GT.llGO  TO  31 

GO  70  3) 

31  G(I1=G« !)-( QI**2/(l.2  7324*AC) ) 
30    AC71=VTS( 11**2/1. 20096E+11 

D  E  LP  5  C  =  RHC  ?  (I )  *  ACT  1  *C  C  RR 
0L7?(  I1=0ELP(I 1-CELPC-CELPEC 
I  FCi.GT.UGO    TO    3J> 
ACTl^GLTPt  I)*Rr'C"1.20096E+U 
AC72  =  li.:,--G(  I  1**2/DI 

33    ACU=DL7PU »*RHCPiI) *1. 200966+11 
ftCT2»12.*G( I1**2/DH 

32  FRIFAGl  I)=ACT1/ACT2 

r 

C    CALCULATE    THE    WILSON    PLOT    PARAMETERS 

AC71=1./RF(  D-.---0.3 

ACT2=1./PR(  11**0.33 

ACT3  =  l./RATr'HU(Il 

X(  I1=ACT1*ACT2*ACT3 
2    CONTINUE 

*E£L*S    TITL  =  <  101/10*'     •/ 

CALL    LSQPL2CM, 1  ,X,  Y  ,  W  I ,  YO  ,0  =  L  Y,  3,  SB  ,  TITLE  ) 
C 
C    DETERMINE    THE    SIEO=R   TATE    CONSTANT 


C 


IF  (N.  EC. 11  GO    TO    96 

CI  =  UN/(PW*2.n*(0H/B<2)l 
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IV  G  LEVEL   21 


MAIN 


DATE 


78C75 


08/46/C9 


GO  TO  95 
9c  CI  =  00/3(2J 
9  5  CONTINUE 

OG  3  1=1, M 

CK.(I)=CI/H20K(I  ) 

3  CONTINUE 

CALCULATE  THE  INSIDE  AND  OUTSIOE  HEAT  TRANSFER  COEFFICIENTS 

DG    4     1-1, M 

IFtN.ECDGG   TO    94 

«C71=CI/DH 

AC72=REl  I)**0.8 

GC    TO    93 
94    ACT1*CI/DI 

4CT2=RE2( i) **0.8 
93    ACT3=PR(I)**( 1./3. ) 

Hid  »*ACTI*ACT2*ACT3*RATMHU(I) 

IFP4.EQ.l)GC    TO    92 

ACTl=P0*3./(AN*UN(I ) ) 

ACT2*(PO/P8AR)*RW 

ACT3  =  f  :/(  PW-HK  Ii  ) 

ACT-V-ACT1-ACT2-ACT3 
HO(I  )=1./ACT4 

GC    TO    91 
92    ACTi  =  CQ/(CI*HI(  I)  I 

ACT2=Y (I)-AW-ACTl 

HC(I  )=i./ACT2 
91    CONTINUE 

CALCULATE    THE    NUSSELT    NUMBER,    THE    STANTCN    NUMBER,     ANC    THE    TUBE    PERFORMANCE 
FACTOR 

IFIN.EC3.DGO    TO    38 

PNU( I )=HI (I )*0H/h20K(  I  ) 

GC  TO  37 
33  PNUII  )=HI (I  )*DI/H20K(  I) 
37  ACT1«PNU(  I  )/?R<  !)**< 1  ./3. ) 

YCCOPC(i)=ACTl/RATMHUm 

ST  (  1  )=P,;  j{  I  )/(?  =  (  I)*PR  (I  )  ) 

C6J(  U=ST(i  )*PR (!)**(  2./ 3.) 

REA'MI)«2.*C3J<  D/FRIFACCIJ 

V  (  l)=V(I)/3600. 

4  CONTINUE 

OUTPUT  SECTION  OF  PROGRAM  STARTS  HERE 

WRI7E(6,204) 
2J4  FORMAT I  '  1'  ,  12X,  'VELOCITY'  ,2  IX,  'UN'  ,24X,'UC  , 24X,'HI  '  ,  24  X  ,  'HO  '  //) 
WRITE  (6,205) 

205  FORMAT  (14X,  '  FT/SEC,  14X,'3TU/(  F*HR*FT**2  )  '  ,  9X,'8TU/(  F*HR*FT**2) 
l',10X,'6TU/(  F*hR*FT**2)  '  ,10X, «STU/<  F*HR*FT**2  )  • /// I 

DC  5  1=1, M 

WRITS (6, 236) VI I >,UN( I ),JC(I) ,HI (I) . HC(I) 

206  F0MA7(9X,F10.2»8X,4(8X,Fl0.3tSX)/J 

5  CONTINUE 
WPITS(6.207 ) 

207  FORMATCO*  ,  3X  ,  •  REYNOLDS    NC'.UX,'     PLAIN    END    REYN    NO', 9X,  'FLOW    RATE 
1    PEA    A?.SA«,7X,'       FRICTION    FACTOR'  ,  10X  ,'  SliDER    TATE    CGNSTANT'//) 

DC    6     1  =  1,  M 

WRITc(6,208)RE< I) ,RE2(I),G(I) ,FRIFAC(IJ ,CK( I ) 

6  CONTINUE 

208  FCRMATtlX,3(3X,F10.2,8X),2(8X,F10.8,dXl/i 
WRITS<6,209> 

209  FORMATCl1  »9X,  "NUSSELT  NO' ,  12X,  •  NU/PR  l/3(  U/UW  )  0  .14'  ,1 IX  ,'  STANTON  N 
IG',17X,'J  FACTOR' , 14X ,'  PERFORM  FACTOR  '//) 

DO  7I»1,M 

WRITE  (6,21 3 )PNU( I) , YCCGRD ( I ) , ST  < I ) , Ce J ( I ) , REANl I ) 

7  CONTINUE 
210FJRMAT{1X,2(8X,F10.5,8X),8X,F10.8,5X,2(3X,F10.7,3X)/) 

WRITE (£,211 ) 
211  FCR'HATi '0'  ,12X,'XIN'  ,13X, 'PRESSURE  DRCP  ( PS  I  )  •  ,  15X,  •  PRANDTL  NO', 19 
IX,  'U/UW'//) 
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"SAN  IV  G  LEVEL   21 


MAIN 


CATE  =  78075 


03/49/09 


>2 
)3 


;7 

:e 

.■9 
.0 

.1 

.2 


00  8  1=1, M 

UaITE(6,212)X(I),0LTP(I),PR(IJ .RATMHUU) 
8  CONTlNuS 
212  FC R v. iT { 5X, E 15. 7, 13X,F 10. 5  ,16X,F10.5,16X,F10.5/) 
DC  20  1=1, M 
ViII*V(I)*3600. 
20  CONTINUE 

IFiN.LT.2JG0  TO  90 
1001  CONTINUE 
STOP 
END 


u  J  J 


APPENDIX  I 

PRESSURE  TAP  TUBE  REDUCTION,  SAMPLE  CALCULATIOI 
FOR  PT-2,  45°  AT  60  PERCENT  FLOW 


In  an  attempt  to  more  accurately  describe  the  actual 
pressure  drop  in  the  three  foot  test  section  of  the  General 
Atomic  tubes,  an  alternate  method  to  that  presented  in  [19]- 
for  predicting  the  pressure  drop  was  sought.   Using  a  specially 
constructed  tube, (see  Figure  14),  several  pressure  drops  were 
measured  across  different  sections  of  the  tube.   These  pres- 
sure drops  are  summarized  in  Tables  21,  22,  and  13  and  were 
used  as  follows  to  arrive  at  a  satisfactory  method  for  pre- 
dicting the  actual  pressure  drop  in  the  three  foot  section  of 
tube. 

To  determine  a  contraction  loss  factor,  the  pressure  drop 
at  the  entrance  must  first  be  determined  using 


AP 


cn 


AP 


6-5 


AP 


Ac 


AP- 


AP 


s  (in  Hg) 


where 


Now 


and 


AP 


AP 


6-5 


Ac 


AP. 


AP 


AP 


AP 


6-5 


Ac 


pressure  drop  measured  between  taps  6  and  5  as 
seen  on  Figure  15 

pressure  drop  due  to   !ie  area  change  at  the  inlet 

pressure  drop  between  tap  5  and  the  beginning 
of  the  test  section 

pressure  drop  in  smooth  section  between  tap  6 
and  the  beginning  of  the  tube. 


13.9  in  Hg 
2 


PV 


TS 


2a 


(  l  -  o 


204 


where 


pV 


TS 


2g, 


and 


Therefore, 

■*   AP 


Ac 


AP. 


10.57  in  Hg, 


'TS 


A 


6.43  in  Hg. 


AP 


3-5 


36  -  L 


3-5 


3-5 


where 


and 


AP 


3-5 


3-5 


Therefore , 

+        AP. 


45.3  in  Hg  , 


31.75  in. 


3.03  in  Hg 


AP 


where 


and 

L 
D 
Therefore, 

-*■         AP 

< 

AP 


cn 


4  fs  G2  (L) 


p2gr  (144) (.4551) 


046 


Re 


0.2 


229 


) 


5.24 


v. 04367 

.048  in  Hg  »  and 

4.42  in  Hg  . 


Knowing  the  inlet  pressure  loss  due  to  contraction  we  can 
determine  the  contraction  loss  factor  from  the  relationship 
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AP 


pV 


en 


en 


K 


TS 


en   2g, 


0.42 


Similarly  at  the  exit,  we  know 


AP 


S  /  • 


APe   =  AP2_6  -  AP5.6  -  4P    -  APT  ♦  APAc-  -ji{1n  Hg) 


where 


AP 


AP 


AP 


2-6 


5-6 


5-3 


AP. 


AP 


Ac 


APs/2 


AP 


60.5  in  Hg  , 

13.9  in  Hg  , 

45.3  in  Hg  , 

3.03  in  Hg  , 

6.43  in  Hg  , 

.024  in  He  ,  and 

4.68  in  Hg 


From 


AP 


pV 


K 


TS 


e      *e  2gc 
it  is  found  that 


K. 


0.44 


The  assumption  that  K  and  K  varied  only  slightly  over  the 
range  of  flows  used  in  this  experiment  was  invoked  as  before 
Thus  the  contraction  loss  factors  and  the  expansion  loss 
factors  were  calculated  at  each  of  the  measured  flow  points, 
and  an  average  of  these  individual  values  was  used  in  deter- 
mining the  non-recoverable  losses  due  to  expansion  and 
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contraction.   The  subscript  2  will  be  used  to  define  the 

friction  factor  and  the  variables  used  to  determine  the 

friction  factor  found  by  using  the  expansion  and  contraction 

loss  factors  as  found  above.   Similarly,  for  the  case  where 

K  „and  K  are  found  using  P9],  the  subscript  1  will  be  used 
en     e 

The  following  constants  are  defined. 


D 


TS 


TS 


D 


h 


3  ft 


0.0241  ft 


lD;TS 


124.5 


'3-5 


3-5 


0.0241  ft 


2.65  ft 


(k) 
^D;3-5 


109.8 


U„  +  Kj,   =   0.33 
en   e'  1 


(Kcn+  Ke]2      =   0.80 


To  determine  the  friction  factor  using  each  of  the 
above  loss  factors,  the  pressure  drop  in  the  test  section 
must  be  found  from: 


where  ; 


(APTS}1,2   "  AP2-6  -  APs  -  (APexp/con>l,2  ^si^ 


AP 


2-6 


27.53  psi  , 


APs    =   .022  psi  , 
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(AP 


exp/con ' 1 


pV 


TS 


29, 


(Kc  +  Ke)1 


1.59  p  s  i  »  and 


(AP 


exp/con'L 


pV 


TS 


2g 


(Kc  +  Ke)2 


3.85  psi  . 


Now  from 


APTS  p  2^c 


and 


3-5 


AP 


3-5 


P  2g, 


4  G'  <*>3-5 


it  can  be  found  that: 


<fTS>l 


0.0108  , 


(fTS)2 


0.0096  ,  and 


(fTS)3.5=   0.0097  . 

It  is  easily  seen  that  the  experimental  results  in- 
dicate that  the  actual  friction  factor  in  the  test  section  is 
less  than  that  computed  using  the  traditional  methods  pre- 
sented in  [19].   It  is  also  noted  that  (^75)3.5  is  very  nearly 

equal  to  (f-rc)?-   Tne  only  rea1  significance  that  can  be 
attached  to  that  fact,  however,  is  that  the  pressure  drop 
was  reduced  to  the  outer  limits  of  the  three  foot  section 
satisfactorily.  Si  nee,  however ,  the  pressure  drop  between  taps 
3  and  5  was  used  in  determining  the  reduction  procedure  - 
the  subsequent  friction  factor  is  coupled  to  fo.c-   There  is 
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obviously  much  more  experimental  work  required  in  this  area 
before  a  universally  accepted  method  can  be  determined. 
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